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1 Introduction 

Novikov's conjecture on the Riemann-Schottky problem: the Jacobians of smooth 
algebraic curves are precisely those indecomposable principally polarized abelian va- 
rieties (ppavs) whose theta-functions provide solutions to the Kadomtsev-Petviashvili 
(KP) equation, was the first evidence of nowadays well-established fact: connections 
between the algebraic geometry and the modern theory of integrable systems is ben- 
eficial for both sides. 

The purpose of this paper is twofold. Our first goal is to present a proof of the 
strongest known characterization of a Jacobian variety in this direction: an indecom- 
posable ppav X is the Jacobian of a curve if and only if its Kummer variety K{X) 
has a trisecant line [36[ 137] . We call this characterization Welters' (trisecant) con- 
jecture after the work of Welters [53]. It was motivated by Novikov's conjecture and 
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Gunning's celebrated theorem [3S]. The approach to its solution, proposed in is 
general enough to be applicable to a variety of Riemann-Schottky-type problems. In 
[24j|38] it was used for a characterization of principally polarized Prym varieties. The 
latter problem is almost as old and famous as the Riemann-Schottky problem but is 
much harder. In some sense the Prym varieties may be geometrically the easiest-to- 
understand ppavs beyond Jacobians, and studying them may be a first step towards 
understanding the geometry of more general abelian varieties as well. 

Our second and primary objective is to take this opportunity to elaborate on 
motivations underlining the proposed solution of the Riemann-Schottky problem, to 
introduce a certain circle of ideas and methods, developed in the theory of soliton 
equations, and to convince the reader that they are algebro-geometric in nature, 
simple and universal enough to be included in the Handbook of moduli. The results 
appeared in this article have already been published elsewhere. 

Riemann-Schottky problem 

Let Hg := {B e Mg{C) = B, lm{B) > 0} be the Siegel upper half space. For 
B eUg let A := As := Z'-i + BZ^ and X := Xb := C^/Ab- Riemann's theta function 

6'(z) := 6'(z, B) ^ e2'^*('"'^)+'^*(™'-^™\ {m, z) = miZi + ■ ■ ■ + rUgZg, (1.1) 

is holomorphic and A-quasiperiodic in z G C^, so 6 := Qb ■= ^^^(0) defines a 
divisor on X. Moreover, (X, [Q]) becomes a ppav, where [9] denotes the algebraic 
equivalence class of 6. Thus Mg/ Sp(2g, Z) ~ Ag, the moduli space of g-dimensional 
ppavs. In what follows we may denote {X, [&]) by X for simplicity. A ppav {X, [&]) G 
Ag is said to be indecomposable if Q is irreducible, or equivalentl}0 if there do not 
exist (Xi, [Qi]) e Ag^ with > 0, i = 1, 2, such that X = Xi x X2 and 9 = 
61 XX2+X1 xGa. 

Let A4g be the moduli space of nonsingular curves of genus g, and let J: Aig 
Ag be the Jacobi map, i.e., for F e Mg, >/(F) is Pic°(F) with canonical polarization 
given by Wg-i = {C e Pic^^^T) \ h"{C) = h^{C) > 0} regarded as a divisor on 
Pic''(F), or more explicitly: taking a symplcctic basis a^, 6^ (z = 1, . . . , g) of i?i(F, Z) 
and a basis cji , . . . , Wg of the space of holomorphic 1-forms on F such that luj — Sij , 
we define the period matrix and the Jacobian variety of F by 

B := (^J^ uj,^ eUg and J(F) := (Xs, [65]) £ , 

respectively. The latter is independent of the choice of (a^, bi). 

J(F) is indecomposable and the Jacobi map J is injective (Torelli's theorem). 
The (Riemann-)Schottky problem is the problem of characterizing the Jacobi locus 
Jg :— J{M.g) or its closure Jg in Ag. For 5 = 2, 3 the dimensions of M.g and Ag 
coincide, and hence Jg = Ag by Torelli's theorem. Since ^4 is of codimension 1 in 
Aa, the case 17 = 4 is the first nontrivial case of the Riemann-Schottky problem. 

A nontrivial relation for the Thetanullwerte of a curve of genus 4 was obtained 
by F. Schottky in 1888, giving a modular form which vanishes on J74, and hence 
at least a local solution of the Riemann-Schottky problem in g = 4, i.e., is an 

^ since principal polarization means parallel translation is the only way to deform 0, translating 
each component of © has the same effect as translating as a whole. 
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irreducible component of the zero locus 1S4 of the Schottky relation. The irreducibility 
of 54 was proved by Igusa [27] in 1981, establishing J4 = ^4, an effective answer to 
the Riemann- Schottky problem in genus 4. 

Generalization of the Schottky relation to a curve of higher genus, the so-called 
Schottky- Jung relations, formulated as a conjecture by Schottky and Jung [54], were 
proved by Farkas- Ranch [3D]. Later, van Geemen [23] proved that the Schottky- Jung 
relations give a local solution of the Riemann-Schottky problem. They do not give 
a global solution when g > 4, since the variety they define has extra components 
already ioi g — 5 (Donagi [TS]V 

More recent development on the Riemann-Schottky problem, as reviewed in [1] 
[51 113). includes a completely new approach of Buser and Sarnak [9 which provides 
an effective way to characterize non-Jacobians. 

Fay's trisecant formula and the KP equation 

Over more than 120 year-long history of the Riemann-Schottky problem, quite a few 
geometric characterizations of the Jacobians have been obtained. Following Mum- 
ford's review with a remark on Fay's trisecant formula [47], and the advent of soliton 
theory and Novikov's conjecture [29l [301 148], much progress was made in the 1980s 
to characterizing Jacobians and Pryms using Fay-like formulas and KP-like equa- 
tions. They are closely related to each other since Fay's formula, written as a biliear 
equation for the Riemann theta function, follows from a difference analogue of the 
bilinear identit'^ 

(f r(<- [fc-i])T(t' + [/c-i])e^(*'-*')'='dfc = 0, (1.2) 

J k—OQ 

which itself is equivalent to the KP hierarchy [lOl E] . Equation (|1.2p can also be 
regarded as a generating function for the Pliicker relations for an infinite dimensional 
Grassmannian. 

Compared with Igusa's work which studies the geometry of 54 and characterize 
the Jacobian locus J^, in this approach Fay-like formulas or KP-like equations are 
used to (in a sense) construct the curve F and thus characterize the Jacobian varieties. 
Therefore this approach to the Riemann-Schottky problem is also related to the 
Torelli theorem; however, the relation is only remote since the conditions like Fay's 
formula and the KP equation contain extra parameters like vector U (and the lack of 
Prym- Torelli does not stop us from studying the Prym-Schottky problem using the 
analogue of this approach) . 

Let us first describe the trisecant formula in geometric terms. The Kummer 
variety K{X) of X G Ag is the image of the Kummer map 

K^Kx: X 3 (e[e, 0](z) | e e ((1/2)Z/Z)f) e CP^'-i (1.3) 

where 0[e,O](z) — 9[e,0]{2z,2B) are the level two theta-functions with half-integer 

^ Here t = {ti,t2, ■ ■ ■) and t' = (t^ , , ■ ■ • ) are two sequences of formal independent variables 
near zero, fc is a formal independent variable near infinity, [k~^] = (1/fc, l/(2fc^), . . . , l/(nfc"), . . . ), 
and T, the so-called tau-function, is a scalar-valued unknown function of the KP hierarchy. For a 
quasiperiodic solution obtained from smooth curve F we have T{t) = e'^('^'>e{J2tiUi + z,_B(F)) for 
some quadratic form Q{t), vectors Ui £ and arbitrary z g C^. Also, Fay's formula itself can in a 
sense be obtained from 111.21 1 by specializing the time variables using the so-called Miwa variables. 
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characteristics e S ((1/2)Z/Z)3, i.e., they equal 9{2{z + Be),2B) up to some expo- 
nential factor so that we have 

e{z + w)9{z - w) = e[e,0](z)e[e,0](w) . (1.4) 

£e((i/2)z/z)s 

We have K{-z) = K{z) and K{X) ~ X/{±1}. 

A trisecant of the Kummer variety is a projective line which meets K{X) at three 
points. Fay's trisecant formula states that if X — J(r), then K(X) has a family of 
trisecants parametrized by 4 points Ai, 1 < i < 4, on F. Namely, identifying a point 
on r with its image under the Abel- Jacobi map F — !> Pic (r) and taking r g Pic"^(r) 
such that 2r ~ A4 — Ai ~ A2 — A3, we have: 

K{r + Ai), K{r + A2) and Kir + A^) are collinear, (1.5) 

i.e., 

^.fA4 + Ai-A2-A3\ ^.f A4- Ai+ A2- A3\ 
2 j' 2 ) 

A are collinear 

if we take the three occurrences of "division by 2" consistent with each other. In 
what follows, the same remark applies if division by 2 in A" appears more than once 
in one formula, as in Theorems 11.251 17.11 

Since we have K{—z) — K{z), condition (|1.5I) is symmetric in all the A^'s. How- 
ever, in its proof as well as its applications the four points tend to play different roles. 
E.g., fixing the 3 points Ai, A2, A3 we may regard it as a one-parameter family of 
trisecants parametrized by A4 or r. Now drop the assumptions that X = J(T) and 
Ai & T C X: suppose A is a ppav such that (|1.5p holds for some Ai, A2, A3 e X 
and infinitely many (hence a one-parameter family of) r £ A. Gunning proved in 
|25| that, under certain nondegeneracy conditions, A is then a Jacobian. 

Gunning's work was extended by Welters who proved that a Jacobian variety can 
be characterized by the existence of a formal one-parameter family of flexes of the 
Kummer variety [63] . A flex of the Kummer variety is a projective line which is 
tangent to K{X) at some point up to order 2. It is a limiting case of trisecants when 
the three intersection points come together. 

In [2] Arbarello and De Concini showed that the assumption in Welters' charac- 
terization is equivalent to a singly infinite sequence of partial differential equations 
contained in the KP hierarchy, and proved that only a first finite number of equa- 
tions in the sequence are sufficient, by giving an explicit bound for the number of 
equations, N = [(3/2)^(7!], based on the degree of K{X). 

Novikov's conjecture 

The second author's answer to Novikov's conjecture [58] illustrated how the soliton 
theory itself can provide natural, useful algebraic tools as well as powerful analytic 
tools to study the Riemann-Schottky problem, as immediately noticed by van der 
Geer [62] , when only an early version of |58j was available: 

An algebraic argument based on earlier results of Burchnall, Chaundy and the 
first author [8l ESI |30] characterizes the Jacobians using a commutative ring R of 
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ordinary differential operators associated to a solution of the KP hierarchy. A simple 
counting argument then shows that only the first 2g + 1 time evolutions in the 
hierarchy are needed to obtain R. Indeed, suppose X = /A appears as an orbit of 
the first 2g + 1 KP flows represented by a "linear motion" (p: C^^"*"^ followed 
by the projection — > X. Then K := keicj) is {g + l)-dimensional, and if (q) e K 
then ^^CidC/ dti — 0, hence by the definition of the KP hierarchy Q — 'Y^^CiPi 
commutes with C. Any two such Q's commute with each other [5S], so the C-algebra 
R' generated by all such Q's is commutative. A simple counting shows that R 
contains an ordinary differential operator of every order n > 2g + 2, which implies 
that R' is maximally commutative and hence R' = R, from the way of constructing 
it. Applying Burchnall et al's theory to R to recover the spectral curve F etc., we 
observe that X ~ J(r). The 2g + 1 KP flows yield a finite number of differential 
equations for the Riemann theta function 9 of X , to characterize a Jacobian. As for 
the number of equations, an easy estimate shows that Ag^ is enough, although more 
careful argument should yield a better bound. Note that this is much smaller than 
Arbarello et al's estimate. 

The analytic tools comes into play when one studies Novikov's conjecture, that 
just the first equation {N = 1!) of the hierarchy, i.e., the KP equation (|1.12l) . 
suffices to characterize the Jacobians: in |58j various tools obtained from analytic 
considerations on the KP equation and family of its solutions were combined with 
the algebraic arguments explained above to prove the conjecture. Even Arbarello and 
De Concini's geometric re-proof of Novikov's conjecture [J used the hardest analytic 
ingredient of [55] as it is, since it had no geometric alternative until Marini's work 
[36] in 1998. Analytic tools are also essential in the proofs of Welters' conjecture and 
its Prym analogue presented in this paper, as condition (C) in each of Theorems ll.61 
[LTQ] [L251 rm Note that (fllO]) . from which condition (C) in Theorem [Tl] follow, 
comes from a generalization of Calogero-Moser system. 

Novikov's conjecture does not give an effective solution of the Ricmann-Schottky 
problem by itself: since it states that X is a Jacobian if and only if 

u = -2{dl In e{Ux + Vy + Wt + Z) + c) 

satisfies ()1.12p for some U, V, W and c, we must eliminate those constants from 
(|1.12|) in order to obtain an effective solution. It is hard to do this explicitly. 

Welters' conjecture 

Novikov's conjecture is equivalent to the statement that the Jacobians are charac- 
terized by the existence of length 3 formal jet of fiexes. In [M] Welters formulated 
the question: if the Kummer variety K{X) has one trisecant, does it follow that X 
is a Jacobian ? In fact, there are three particular cases of the Welters conjecture, 
corresponding to three possible configurations of the intersection points (a, &, c) of 
K{X) and the trisecant: 

(i) all three points coincide {a — b — c); 

(ii) two of them coincide {a — b ^ c); 

(iii) all three intersection points are distinct {a b =^ c ^ a). 
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Of course the first two cases can be regarded as degenerations of the general case (iii). 
However, when the presense of only one trisecant is assumed, all three cases are inde- 
pendent and require separate treatment. The proof of case (i) of Welters' conjecture 
was obtained by the first author in |36j : 

Theorem 1.6 An indecomposable principally polarized abelian variety {X,9) is the 
Jacobian variety of a smooth algebraic curve of genus g if and only if there exist 
g-dimensional vectors U ^ 0,V,A , and constants p and E such that one of the 
following three equivalent conditions are satisfied: 
{A) the equality 

{dy-dl+u)i^^Q, (1.7) 

where 

holds, for an arbitrary vector Z ; 

(B) for all theta characteristics e E (^Z/Z)^ 

{dv ~dlj~2pdu + {E- p2)) e[e, Q]{A/2) = 

(here and below djj , dv are the derivatives along the vectors U and V , respectively). 

(C) on the theta-divisor Q = {Z e X \ 0{Z) = 0} 

[{dvOf - {d^e)^]d^e + 2[dledle - dvedudv0]du9 + [die - d^e] {duOf = (mod e) 

(1.9) 

The direct substitution of the expression (|1.8p in equation (|1.7I) and the use of the 
addition formula for the Riemann theta-functions shows the equivalence of conditions 
{A) and {B) in the theorem. Condition {B) means that the image of the point A/2 
under the Kummer map is an inflection point (case (i) of Welters' conjecture). 

Condition (C) is the relation that is really used in the proof of the theorem. 
Formally it is weaker than the other two conditions because its derivation does not 
use an explicit form (|1.8p of the solution -0 of equation (|1.7I) . but requires only 
an existence of a meromorphic solution: consider a holomorphic function t{x, y) 
of a complex variable x depending smoothly on a parameter y, and assume that 
in a neighborhood of a simple zero ry(?/) of function r (that is, T{r]{y),y) = and 
dxT{r](y),y) ^ 0) equation (|1.7p with potential u = —2d^\nT has a meromorphic 
solution -0. Then the equation 

ii^2w, (1.10) 

holds, where the "dots" denote derivatives in y, and w is the third coefficient of the 
Laurent expansion of the function u at the point 77, i.e., 

2 

u{x, y) = ---2 + v{y) + w{y){x - i]{y)) + ■■■ . 

Equations (|1.10l) was first derived in H| where the assertion of the theorem was proved 
under the assumptior|f| that the closure of the group in X generated by A coincides 



^under different additional assumptions the corresponding statement was proved in the earUer 
works [33ll46l 
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with X. Expanding the function in a neighborhood of a point z £ 8 := {z | 0{z) = 
0} such that du9{z) ^ 0, and noting that the latter condition holds on a dense subset 
of 6 since B is indecomposable, it is easy to see that equation (ll.lOp is equivalent to 

dm). 

Equation (|1.7p is one of the two auxiliary linear problems for the KP equation. 
Namely, the compatibility condition of (|1.7p and the second auxiliary linear equation 

dt-dl + ^ud^ + u^^^Q (1.11) 



is equivalent to the KP equation [IHl IHi] : 



3 d f I 3 , 

-^Uyy = —\Ut - ^U^^XX - -jUU^ ] . (1-12) 



For the first author, the motivation to consider not the whole KP equation but just 
one of its auxiliary linear problem was his earlier work 33^ on the elliptic Calogero- 
Moser (CM) system, where it was observed for the first time that equation (|1.7p is 
all what one needs to construct the elliptic solutions of the KP equation. Moreover, 
the construction of the Lax representation with a spectral parameter and the corre- 
sponding spectral curves of the elliptic CM system proposed in [33] can be regarded 
as an effective solution of the inverse problem: how to reconstruct the algebraic curve 
from the matrix B if its Kummer variety admits one fiex with the vector U (in the 
assumption of the Theorem) which spans an elliptic curve in the abelian variety X . 
Briefiy, that solution of the reconstruction problem can be presented as follows: 
// the vector U spans an elliptic curve E C X , then the equation 

0{Ux + Vy + Z) = (1.13) 

for a generic Z has g simple roots Xi{y) depending on y (they are just intersection 
points of the shifted elliptic curve E + Vy + Z C X with the theta-divisor O (Z X). 
These roots define g x g matrix L{y^ z) with entries given by 

Lu{t,z)^^Xi, ^ ^{x^ - Xj,z), i^j, (1-14) 



where 



^x,z):^ ^1^, ^'l e^^^-^^, (1.15) 



(j(z)a{x) 

with C and a the standard Weierstrass functions. 

The spectral curve Pcm of the CM system is the normalization at the point k = 
oo, z = of the closure in x E of the affine curve given in C x (E \0) by the 
characteristic equation 

R{k,z)^detikl + L{y,z)) = 0. (1.16) 

Under the assumptions of the theorem, the CM curve Pcm does not depend on y and 
is the solution of the inverse problem. 

Without an assumption on U the proof of Theorem 11.61 is much more complex 
and less effective. The ultimate goal is to construct, under the assumption that 
the condition (C) is satisfied, a ring of commuting ordinary differential operators. 
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because, as shown in [S], a pair of commuting differential operators Li, L2 satisfies an 
algebraic relation R{Li,L2) = 0. This is the key moment, when an algebraic curve 
emerges in the proof. It then remains only to show that the corresponding curve is 
the solution of the inverse problem. 

The first step in the proof is to introduce in the problem a formal spectral param- 
eter. It is analogous to the introduction of the spectral parameter in the Lax matrix 
for the elliptic CM system. This parameter k appears in the notion of a formal wave 
solution of equation (|1.7I) . 

The wave solution of (11.71) is a solution of the form 

^{x, y, k) = efc-+(*^'+'')!^ (1 + £ y) ■ (1-17) 

The aim is to show that under the assumptions of the theorem there exists a unique, 
up to multiplication by a constant factor c(fc), formal wave solution such that 

^ _ Ts{Ux + Vy + Z,y) 



6{Ux + Vy + Z) 



where Tg{Z,y), is an entire function of Z. 

As it was stressed above, strictly speaking the KP equation and the KP hierarchy 
are not present in the assumptions of the theorem, but the analytical difficulties in 
the construction of the formal wave solutions of (|1.7p can be traced back to those in 
the second author's proof [58] of Novikov's conjecture. 

The main idea of proof in [55] is to show that if tq = e'^^ ^^9{Ux + Vy + Wt + Z) 
satisfies the KP equation in Hirota's for nfl 

{Dl + iDl~^D,Dt)To-To^Q, 

so that u — — 29^To satisfies the KP equation (|1.12l) . then it can be extended to 
a T-function of the KP hierarchy, as a global holomorphic function of the infinite 
number of variables t — (ti) — (ti, ^2, ^3, • • • )j with ti = x, t2 = y, ~ t. Local 
existence of r directly follows from the KP equation. The global existence of the 
r-function is crucial. The rest is a corollary of the KP theory and the theory of 
commuting ordinary differential operators developed by Burchnall-Chaundy [5] and 
the first author [^150]. 

The core of the problem is that there is a homological obstruction for the global 
existence of r. It is controlled by the cohomology group if^(C^\E, V), where singular 
locus E is defined as 9t/-invariant subset of the theta-divisor 9 and V is the sheaf of 
^(/-invariant meromorphic functions on \ S with poles along Q. The hardest part 
of [SS] , as clarified in [3] , is the proof that the locus S is empty 0. 

The coefficients of the wave function are defined recurrently by the equation 
2du£,s+i — dy£^s — dijS^s + u^s- It turned out that equation (|1.9|) in the condition (C) 
of the theorem are necessary and sufficient for the local existence of meromorphic 
solutions. The global existence of is controlled by the same cohomology group 
H^{C^ \ S, V) as above. Fortunately, in the framework of our approach there is 



* We define P(D^, ...)/■/:= P{d^>,. . .){f{x + x' , . . :)f{x - x', . . . ))|^/ = ...=o for a polynomial 
or a power series P; a Hirota equation is an equation of the form P{Dx, ...)/■/ = 0; see 1111 158| . 

^The first author is grateful to Enrico Arbarello for an explanation of these deep ideas and a 
crucial role of the singular locus S, which helped him to focus on the heart of the problem. 
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no need to prove directly that the bad locus is empty. The first step is to construct 
certain wave solutions outside the bad locus. We call them A-periodic wave solutions. 
They are defined uniquely up to ^[/-invariant factor. The next step is to show 
that for each Z ^ Yi the A-periodic wave solution is a common eigenfunction of 
a commutative ring of ordinary difference operators. The coefficients of these 
operators are independent of ambiguities in the construction of ?/'• For the generic 
Z the ring is maximal and the corresponding spectral curve F is independent. 
The correspondence j : Z i — > A^ and the results of the works [H [29l |30l |48] , where 
a theory of rank 1 commutative rings of differential operators was developed, allows 
us to make the next crucial step and prove the global existence of the wave function. 
Namely, on {X \ S) the wave function can be globally defined as the preimage j*ipBA 
under j of the Baker- Akhiezer function on F and then can be extended on X by 
usual Hartogs' arguments. The global existence of the wave function implies that 
X contains an orbit of the KP hierarchy, as an abelian subvariety. The orbit is 
isomorphic to the generalized Jacobian J(F) = Pic°(r) of the spectral curve ([58]). 
Therefore, the generalized Jacobian is compact. The compactness of J(F) implies 
that the spectral curve is smooth and the correspondence j extends by linearity and 
defines the isomorphism j : X ~> ^(F). 

The proof of Welters' conjecture was completed in [37]. First, here is the theorem 
which treats case (ii) of the conjecture: 

Theorem 1.19 An indecomposable, principally polarized abelian variety {X,9) is 
the Jacobian of a smooth curve of genus g if and only if there exist non-zero g- 
dimensional vectors U ^ A (mod A), V, such that one of the following equivalent 
conditions holds: 

(A) The differential-difference equation 

{dt -T + u{x,t))ip{x,t) ^0, r = e^- (1.20) 

is satisfied for 

u= {T -l)v{x,t), V ~dtln0{xU + tV + Z) (1.21) 

and 

^^0{A^xU + tV + Z) 

^ 0{xU + tV + Z) ' ^ ' 

where p, E are constants and Z is arbitrary. 

(B) The equations 

dvQle, 0] {{A - U)/2) - ef e[e, 0] {{A + U)/2) + EQ[e, 0] {{A - U)/2) = 0, 

are satisfied for all e £ (^Z/Z)^. Here and below dv is the constant vector field on 
corresponding to the vector V . 

(C) The equation 

dv [0{z + u) e{z - u)] dvO{z) = [e{z + u) e{z - u)] dlyeiz) (mod e) (1.23) 

is valid on the theta-divisor Q = {Z G X \ 9{Z) = 0}. 
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Equation (|1.20l) is one of the two auxiliary linear problems for the 2D Toda lattice 
equation 

d^d^^^ = e'^-.-i-'^. _ e'^.-'^.+i^ (1.24) 

which can be regarded as a partial discretization of the KP equation. The idea to 
use it for the characterization of the Jacobians was motivated by |36| and the first 
author's earlier work with Zabrodin |45) . where a connection of the theory of elliptic 
solutions of the 2D Toda lattice equations and the theory of the elliptic Ruijsenaars- 
Schneider system was established. In fact, Theorem 1 1.1 91 in a slightly different form 
was proved in [45| under the additional assumption that the vector U spans an elliptic 
curve in X. 

The equivalence of {A) and [B) is a direct corollary of the addition formula for 
the theta-function. The statement {B) is the second particular case of the trisecant 
conjecture: the line in CP^ ~^ passing through the points K{{A — U)/2) and K{{A + 
U)/2) of the Kummer variety is tangent to K{X) at the point K{{A — U)/2). 

The affirmative answer to the third particular case, (iii), of Welters' conjecture is 
given by the following statement. 

Theorem 1.25 An indecomposable, principally polarized abelian variety {X,9) is the 
Jacobian of a smooth curve of genus g if and only if there exist non-zero g-dimensional 
vectors U^V=^A=/=U( mod A) such that one of the following equivalent conditions 
holds: 

(A) The difference equation 

ip{m, n + 1) = i/;(m + 1, n) + u{m, n)^{m, n) (1-26) 

is satisfied for 

_ e{{m + 1)U + {n + 1)V + Z) ejmU + nV + Z) 

~ e{mU + {n+l)V + Z)9{{m + l)U + nV + Z) ^' ' 

and 

0[mU + nV + Z) 
where p, E are constants and Z is arbitrary. 

(B) The equations 

A^U~V\ , ^,/'A + U~V\ ^.fA + V-U 



e[e,0] j +e^e[e,0] ( ) =e-e[e,0] 

are satisfied for all e G (^Z/Z)^. 
(C) The equation 

0{z+u)e{z-v)e{z-u+v)+0{z-u)e{z+v)9{z+u-v) = (mode) (1.29) 

is valid on the theta-divisor O = {Z G X \ 9{Z) = 0}. 

Under the assumption that the vector U spans an elliptic curve in X, Theorem ll.251 
was proved in , where the connection of the elliptic solutions of BDHE and, the 
so-called, elliptic nested Bethe Ansatz equations was established. 
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Equation (|1.26l) is one of the two auxiliary linear problems for the so-called bilinear 
discrete Hirota equation (BDHE): 

Tn{l + 1, m)T„(/, m + 1) - r„(/, m)T„{l + 1, m + 1) + r„+i(/ + 1, m)T„-i{l,m + 1) = 

(1.30) 

At the first glance all three nonlinear equation: the KP equation, the 2D Toda 
equation, and the BDHE equation, look quite unlikely. But in the theory of integrable 
systems it is well-known that these fundamental soliton equations are in intimate 
relation, similar to that between all three cases of the trisecant conjecture. Namely, 
the KP equation is as a continuous limit of the BDHE, and the 2D Toda equation 
can be obtained in an intermediate step. 

The structure of the statements of the last two theorems, and the structure of 
their proofs look almost literally identical to that in Theorem 11.61 To some extend 
that is correct: in all cases the first step is to construct the corresponding wave 
solution. The conditions (C) in all three cases play the same role. They ensure the 
local existence of the wave function. The key distinction between the differential and 
the difference cases arises at the next step. As it was mentioned above, in the case 
of differential equations a cohomological argument [SHI Lemma 12] can be applied 
to glue local solutions into a global one. In the difference case there is no analog of 
the cohomological argument and we use a different approach. Instead of proving the 
global existence of solutions we, to some extend, construct them by defining first their 
residue on the theta-divisor. It turns out that the residue is regular on Q outside the 
singular locus E. Surprisingly, it turns out that in the fully discrete case the proof 
of the statement that the singular locus is in fact empty can be obtained at much 
earlier stage than in the continuous or semi-continuous case. In part, it is due the 
drastic simplification in the fully discrete case of the corresponding equation on the 
theta-divisor (compare (|1.29l) with (11.91) ). 

Structure of the article 

In the next section we introduce the basic concept of the algebro-geometric integra- 
tion theory of soliton equation, that is the concept of the Baker- Akhiezer function, 
which is defined by its analytic properties on an algebraic curve with fixed local co- 
ordinates at marked points. The uniqueness of the Baker- Akhiezer function implies 
that it is a solution of certain linear differential equations. The existence of the 
Baker- Akhiezer function is proved by explicit theta-functional formula, which then 
leads to explicit theta-functional formulae for the coefficients of the corresponding 
equations. That proves "the only if" part in all the theorems above. 

In section 3, we introduce the KP hierarchy in Sato's form as a system of com- 
muting flows on the space of formal pseudodifferential operators. Because, the flows 
commute, the hierarchy can be reduced to the stationary points of one of the flows 
(or their linear combination). That is a reduction from a spatially two-dimensional 
system to a spatially one-dimensional systemic] Under this reduction, the KP hier- 

® Here the term "spatially two-dimensional (resp. one-dimensional) system," also known as "sub- 
subholonomic (resp. subholonomic) system" or "(2 -|- l)-d (resp. (1 -|- l)-d) system," means the one 
whose "general solution" depends on functions of two variables (resp. one variable), or equivalently, 
on doubly infinite (resp. singly infinite) sequences of parameters in the formal power series set-up. 
(The word "space" is associated to the notion of free parameters because in an initial value problem 
of a partial differential equation the free parameters for a solution are given by its initial data, which 
are given on a "space-like" hypersurface.) E.g., since initial data for the KP hierarchy, i.e., C\t=o 



11 



archy defined first on "a space" of infinite number of functions of one variable (tlie 
coefficients of a pseudodifferential operator) is equivalent to a system of commuting 
flows on the space of finite number of functions of one variable. For the case of 
stationary points of a linear combination of the first n flows of the KP hierarchy 
these functions are coefficients of a differential operator Ln of order n. One may take 
one step further and consider stationary points of two commuting flows. It turns 
out that if the corresponding integers n and m are co-prime, then the corresponding 
orbits of the whole hierarchy are finite-dimensional and can be identified with certain 
subspaces of the finite-dimensional linear space of solutions to the system of ordinary 
difTerential equations: 

n— 1 m — 1 

[L„, Lm] =0, L„ = 9: + J2 u^{x)^l, Lm = + v,{x)di (1.31) 

i=0 j=0 

This is a setup explaining the role of commuting operators in the modern theory of 
integrable systems. 

As a purely algebraic problem it was considered and partly solved in the remark- 
able works of Burchnall and Chaundy [8 in the 1920s. They proved that for any pair 
of such operators there exists a polynomial in two variables such that R{L„^ Lm) = 0. 
Moreover, they proved that if the orders n and m of these operators are co-prime, 
(n^m) — 1, and the algebraic curve F defined in by equation R(X, ^) = is 
smooth, then the commuting operators are uniquely defined by the curve and a set 
of g points on F, where g is the genus of F. In such a form, the solution of the 
problem is one of pure classification: one set is equivalent to the other. Even the 
attempt to obtain exact formulae for the coefficients of commuting operators had not 
been made. Baker proposed making the programme effective by looking at analytic 
properties of the eigenfunction ^. The Baker program was rejected by the authors 
of [5] consciously (see the postscript of Baker's paper [5]) and all these results were 
forgotten for a long time. 

The theory of commuting differential operators and its extension to the difference 
case is presented in Section 4. The outline of the proof of the trisecant conjecture 
is in Section 5. In Section 6 we present a solution of the characterization problem 
for Prym varieties which was obtained by Grushevsky and the first author f [381124) 1. 
The last Section 7 is devoted to a theory of abelian solutions of the soliton equation. 
The notion of such solutions was introduced by the authors in [42l |43] , where it was 
shown that all of them are algebro-geometric. The theory of abelian solutions can be 
regarded as an extension of the results above to the case of non-principally polarized 
abelian varieties. 

2 The Baker- Akhiezer functions — General scheme 

Let F be a nonsingular algebraic curve of genus g with N marked points Pa and 
fixed local parameters k~^{Q) in neighborhoods of the marked points. The basic 

for C in 114.411 . are given by a singly infinite sequence of one-variable functions {vs{x)}g—i^2,... or, by 
expanding each Va{x) in a power series Vs{x) = doubly infinite sequence of parameters 

{Dgjjs—i 2,...;i=o,l,...! the KP hierarchy is a "spatially two-dimensional system." For 2 < n S Z the 
n-reduction of the KP hierarchy (KdV if n = 2, Boussinesq if n = 3, etc.) is defined by imposing 
the condition that £" is a differential operator. Since, as an ordinary differential operator, C"\t—o 
depends on finite number of one-variable functions and hence on finite number of singly-infinite 
sequences, it is a "spatially one-dimensional system." 
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scalar multi-point and multi-variable Baker-Akhiezer function ^jJ(t, Q) is a function 
of external parameters 

t = {tea), a = 1, . . . , A^; i = 0, . . . ; ^ to.,o = 0, (2.1) 

a 

only finite number of which is non-zero, and a point Q £ T . For each set of the 
external parameters t it is defined by its analytic properties on T. 

Remark. For the simplicity we will begin with the assumption that the variables 
ta_n are integers, i.e., ta_n G 

Lemma 2.2 For any set of g points 71, . . . ,7g in a general position there exists a 
unique (up to constant factor c(t)) function ip(t,Q), such that: 

(i) the function ip (as a function of the variable Q GT) is meromorphic everywhere 
except for the points and has at most simple poles at the points 71 , . . . , 7^ ( if all 
of them are distinct); 

(ii) in a neighborhood of the point the function ip has the form 



(2.3) 



where k^ — ka{Q) is the reciprocal of a local parameter at Pa, i.e., k^^ G mp^ \™p„ • 

From the uniqueness of the Baker-Akhiezer function it follows that: 

Theorem 2.4 For each pair (a, n > 0) there exists a unique operator La^n of the 
form 

n-l 

i«.n = a2,i + E4"'"^W<i' (2-5) 

(where da^n ~ d/dta^n) such that 

{da,n-La,n)Ht,Q)=0. (2.6) 

The idea of the proof of the theorems of this type proposed in [29], [30] is universal. 

For any formal series of the form (12. 3p their exists a unique operator La,n of the 
form (|2.5p such that 

{da,n - La,n) Q) = 0{k-^) CXp ta,^kC\ ■ (2.7) 

The coefficients of L^^n are universal differential polynomials with respect to (,s,a- 
They can be found after substitution of the series (|2.3p into (|2.7p . 

It turns out that if the series (12. 3|) is not formal but is an expansion of the 
Baker-Akhiezer function in the neighborhood of Pa the congruence (|2.7p becomes an 
equality. Indeed, let us consider the function ipi 

^l = {da,n-La,n)^Ai,Q). (2.8) 

It has the same analytic properties as V' except for the only one. The expansion of 
this function in the neighborhood of Pa starts from 0(k~^). From the uniqueness of 
the Baker-Akhiezer function it follows that -01 = and the equality (12. 6p is proved. 
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Corollary 2.9 The operators La^n satisfy the compatibility conditions 

Remark. The equations p.lOp are gauge invariant. For any function c{t) operators 

La,n = cLa^nC ^ + {da.nC)c ^ (2.11) 

have the same form (|2.5p and satisfy the same operator equations (I2.10p . The gauge 
transformation (j2.1ip corresponds to the gauge transformation of the Baker- Akhiezer 
function 

iP{t,Q) = c{t)^{t,Q) (2.12) 



In addition to differential equations (12.61) the Baker-Akhiezer function satisfies an 
infinite system of differential-difference equations. Recall that the discrete variables 
ta,Q are subject to the constraint J2a ^".o — 0- Therefore, only the first {N — 1) 
of them are independent and t^r^o = ~S^=i^^a,o- Let us denote by Tq, a — 
1, . . . , TV— 1, the operator that shifts the arguments t^.o — ^ ia,o+l and tN^ tj^fl — ^, 
respectively. For the sake of brevity in the formulation of the next theorem we 
introduce the operator T/v = T^^. 

Theorem 2.13 For each pair (a, n > 0) there exists a unique operator L^^n of the 
form 

n-l 

= + .r'"^(t) = 0. (2.14) 

3=0 

such that 



£a.n) ^(t,g) =0. (2.15) 



The proof if identical to that in the differential case. The operators La^n are defined 
by congruence insuring that the resulting function satisfies all the condition of the 
Baker-Akhiezer function plus vanishing of one of the leading coefficients. After that 
the uniqueness of the Baker-Akhiezer function implies that the congruence is in fact 
the equality. 

Corollary 2.16 The operators L^.n satisfy the compatibility conditions 

It should be emphasized that the algebro-geometric construction is not a sort of 
abstract "existence" and "uniqueness" theorems. It provides the explicit formulae 
for solutions in terms of the Riemann theta-functions. They are the corollary of the 
explicit formula for the Baker-Akhiezer function: 

Theorem 2.18 The Baker-Akhiezer function is given by the formula 

V'(i, P) = C{t) exp (^2^ ^a,^^^a,^(-P) j e{A{P) + Z) ' ^ ^ ^ 

Here the sum is taken over all the indices (a, i > 0) and over the indices (a, 0) with 
a — \, . . . ^ N — 1, and: 
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a) ^a,iiP) is the abelian integral, fta^i(P) = J dflaj, corresponding to the unique 
normalized, dfla,i = 0, meromorphic differential on T , which for i > has the 
only pole of the form dila,i — d [k]^ + 0(1)) at the marked point Pa and for i = Q 
has simple poles at the marked point Pa and P^ with residues ±1, respectively; 

h) 2'KiUa,j is the vector of b-periods of the differential dflaj, i-C, 

c) A{P) is the Abel transform, i.e., a vector with the coordinates A{P) = cLok 

d) Z is an arbitrary vector (it corresponds to the divisor of poles of Baker- Akhiezer 
function). 

Notice, that from the bihnear Riemann relations it fohows that the expansion of the 
Abel transform near the marked point has the form 

oo ^ 

AiP)^AiPa)-y^-Ua.^ka' (2.20) 



Example 1. One-point Baker- Akhiezer function. KP hierarchy 

In the one-point case the Baker- Akhiezer function has an exponential singularity at 
a single point Pi and depends on a single set of variables ti = ti^i. Note that in this 
case there is no discrete variable, ti^Q = 0. Let us choose the normalization of the 
Baker- Akhiezer function with the help of the condition o = 1, i.e., an expansion of 
tjj in the neighborhood of Pi equals 



^'(ii, . . . , Q) = exp ^ t,k^ ^ + E ^'^(^)^"' • (2.21) 
Under this normalization (gauge) the corresponding operator L„ has the form 

n-2 

Ln^d'^ + Y.^^^dl. (2.22) 



i=0 

32 



For example, for n — 2,3 after redefinition x ~ ti we have L2 — ~ u,L3 ~ 
— ^udx — w with 

u{x,t2,...) = 2dx^i{x,t2,...), (2.23) 



Therefore, if we define y — t2,t — t^, then u{x, y, t, t^, . . .) satisfies the KP equation 

The normalization of the leading coefficient in (|2.2ip defines the the function 
c{t) in (|2.19p . That gives the following formula for the normalized one-point Baker- 
Akhiezer function: 

M Q) ~ exp (Tt n (P)] mP)+j:U.U + Z)9iZ) 

(shifting Z if needed we may assumed that ^(Pi) = 0). In order to get the explicit 
theta-functional form of the solution of the KP equation it is enough to take the 
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derivative of the first coefficient of the expansion at the marked point of the ratio of 
theta- functions in the formula (|2.24[) . 

Using (I2.20[) we get the final formula for the algebro-geometric solutions of the 
KP hierarchy |30] 

oo 

u{ti,t2,...) = ~2df\n 9 C^U^U + Z) + const. (2.25) 



Example 2. Two-point Baker- Akhiezer function. 2D Toda hierarchy 

In the two-point case the Baker-Akhiezer function has exponential singularities at 
two points Pa, a — 1,2, and depends on two sets of continuous variables ta.i>o- In 
addition it depends on one discrete variable n = ti^ = — t2,o- Let us choose the 
normalization of the Baker-Akhiezer function with the help of the condition ^i^o = 1, 
i.e., in the neighborhood of Pi the Baker-Akhiezer function has the form: 

^{n, to..>o, Q) = fc" exJf^ ^i^'^O (l + £ ^)'^r') ' (2.26) 

and in the neighborhood of P2 

^{n, ta^^>Q,Q) = ^2"" cxp (f^ h,^k{) (f2 6,s (", ^^1"') , (2-27) 



s=0 



According to Theorem l2.4[ the function ■0 satisfies two sets of differential equations. 
The compatibility conditions (|2.10l) within the each set can be regarded as two copies 
of the KP hierarchies. In addition the two-point Baker-Akhiezer function satisfies 
differential difference equation (|2.14p . The first two of them have the form 

(5i,i -r + w)V^ = 0, (92,1 - wT-i)^/> = 0, (2.28) 

where 

M= (r-l)ei,i(n,t), u; = e^"-^"-\ e-^-^*) =e2,o(ri,t) (2.29) 

The compatibility condition of these equations is equivalent to the 2D Toda equation 
(|1.24p with ^ = ^1^1 and r] = ^2,1- The explicit formula for 0„ is a direct corollary of 
the explicit formula for the Baker-Akhiezer function. The normalization of ip as in 
(pr^ defines the coefficient c in (PTTOl) 



e{A{P) +nU + Y: Ua,^ta,^ + Z) 9{Z) 

0{nU + Y: UaJo.,, + Z) {0{A{P) + Z) ■ 



1p = exp Inilifi + ta,iila,i{P) ) i 4 i r^\fQ^ A^u\ i ' (2-30) 



If we denote a; = 0, t — ti.i and set ti,i>i = i2,i>o = Oj then up to a constant in 
{x,t) factor the formula (j2.30p coincides with (|1.22p . Expanding ip at Pi we get the 
formula for the coefficient u in in the first linear equation ()2.29p . which coincides 
with (|1.2ip . That proved "the only if" part of Theorem 11.191 

Example 3. Three-point Baker-Akhiezer function 

Starting with three-point case, in which the number of discrete variables is 2, the 
Baker-Akhiezer function satisfies certain linear difference equations (in addition to 
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the differential and the differential-difference equations (I2.6p . (l2.15p V The origin of 
these equations is easy to explain. Indeed, if all the continuous variables vanish, 
ta,i>o — 0, then the Baker- Akhiezer function i^n.miP), where n = —ti^^ m = —^2,01 
is a meromorphic function having pole of order n + to at P3 and zeros of order n and 
m at Pi and P2 respectively, i.e., 

^n,,n(^H^{D + n{P3-Pi)+m{P3-P2)), i? = 7i + .-.+73 (2.31) 

The functions ijjn+i,m,'4'n,m+i,'4'n,m are all in the linear space H^(D -I- (n -I- m -I- 
1)^3 — nPi — mP2). By Riemann-Roch theorem for a generic D the latter space is 2- 
dimensional. Hence, these functions are linear dependent, and they can be normalized 
such the the linear dependence takes the form (|1.26p . The theta- functional formula 
for the Baker- Akhiezer function directly implies formulae (|1.27p . (|1.28l) and proves 
"the only if" part of Theorem [OSl 

For the first glance it seems that everything here is within the framework of 
classical algebraic-geometry. What might be new brought to this subject by the 
soliton theory is understanding that the discrete variables ta^Q can be replaced by 
continuous ones. Of course, if in the formula (j2.19p the variable tafi is not an integer, 
then ip is not a single valued function on F. Nevertheless, because the monodromy 
properties of ip do not change if the shift of the argument is integer, it satisfied the 
same type of linear equations with coefficients given by the same type of formulae. 
It is necessary to emphasize that in such a form the difference equation becomes 
functional equation. 

Remark. In the four-point case there is three discrete variables n, m, I. In each two 
of them the Baker- Akhiezer function satisfies a difference equation. Compatibility of 
these equations is the BDHE equation (|1.30p . 



3 Dual Baker- Akhiezer function 

The concept of the dual Baker- Akhiezer function ip^{t,P) is universal and is at the 
heart of Hirota's bilinear form of soliton equations, and plays an essential role in our 
proof of Welters' conjecture. It is necessary to emphasize that, although the concept 
is universal, the definition of the dual Baker- Akhiezer function depends on a choice 
of dual divisor — it + ■ ■ ■ + 1^ ■ As it will be shown later the notion of duality 
between divisors of t/i and ■0^ reflects a choice of one of the variables t^fi or ta,i- In 
all the cases the pole divisor Z?+ of the dual Baker-Akhiezer function is defined by 
the equation 

D + D+ ^K + n(^J{T) (3.1) 

where K \s a. canonical class and k is a certain degree 2 divisor, that encodes the 
type of duality. Depending on its choice, the dual Baker-Akhiezer function is then 
defined by the following analytic properties: 

i) the function ■0"'" (as a function of the variable P € F) is meromorphic everywhere 
except for the points Pa and has at most simple poles at the points j • ■ • : 7^ (if 
all of them are distinct); 

(ii) in a neighborhood of the point Pa the function ?/' has the form 

^(t,Q) = /c-*-°exp[f;-i„,,fcj,Vf;^+,(t)/c-^Y ka^kaiQ). (3.2) 



17 



In fact it is the same Baker-Akhiezer type function and, therefore, admits the same 
type of exphcit theta-function formula: 

The basic type of duality and their meaning are explained below in two examples. 

Example 1. One-point case. Duality for a continuous variable 

The notion of dual Baker-Akhiezer function in the one point case was first introduced 
in [TU]. In this dual divisor is defined by p.ip where k = 2Pi. In other words, for 
a generic effective degree g divisor D there exists a unique meromorphic differential 
dil with pole of degree 2 at Pi, dil = d{ki + 0(1)) having zeros at the points 7^; in 
addition it has g more zeros that are denoted by 7^^ , ■ • ■ , 7^ • 

The functions "0 and ip~^s have essential singularities, their product or products 
of their derivatives are meromorphic functions on F. Moreover, from the definition 
of the duality it follows that after multiplication by corresponding differential dfi 
one gets a meromorphic differential on F with the only pole at Pi. That proves the 
following statement. 

Lemma 3.4 Let ip and ip'^ be the Baker-Akhiezer function and its dual. Then the 
following equations hold: 

Tesp,{iP+idi^))dn^O, J =0,1,.... (3.5) 

Equations (13. Sp allows to express the coefficients ^+ of the expansion of the dual 
function ip^ at Pi as universal differential polynomials in terms of the coefficients 
^s' of the Baker-Akhiezer function. The first such equation is Ci + 'fi^ = 0- Another 
corollary of p.Sp is infinite number of bilinear identities for the theta-function, that 
one obtains after substitution of (I2.19p . (|3.3p into (|3.5p . These identities are usually 
called Hirota's bilinear equations. 

Corollary 3.6 Let be the Baker-Akhiezer function and Li be the linear operator 
of the form I12.22\) such that (On — Ln)ip — 0. Then the dual Baker-Akhiezer function 
is a solution of the formal adjoint equation 

4>+id„-L„)^0 (3.7) 

Recall that the right action of a differential operator is defined as a formal adjoint 
action, i.e., f^di — ~dif^ (and the left-hand side of this formula should not be 
confused with the more common differentiation-followed-by-multiplication construc- 
tion for a differential operator). The proof of the corollary will be given in the next 
section. 

Example 2. Two-point case. Duality for a discrete variable 

In the two-point case, in which there is one discrete variable n, the dual divisor D'^ 
is defined by p.ip with k = Pi -I- P2, i.e., 7s and 7^' are zeros of a differential dVL 
having simple poles at the marked points Pi and P2. Without loss of generality we 
may assume that at these points it has residues =Fl- 
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Lemma 3.8 Let tp and be the Baker- Akhiezer function and its dual. Then the 
following equations hold: 

respi (V'+(T») df7 = 0, i = l,2.... (3.9) 

By definition of the duality, the differential on the left-hand side of p.9p has pole 
only at Pi. Hence its residue vanishes. Note also that the differential ip^ipdn has 
poles at Pi and P2 . The constant c"*" in the normalization of the dual Baker- Akhiezer 
function is chosen such that 

respi (ip+ip) dn = 1. (3.10) 



Corollary 3.11 Let ^ be the Baker- Akhiezer function and let Ln be the linear op- 
erator of the form 

n-l 

= r + (3.12) 

such that (di^i — Ln)^} = 0. Then the dual Baker- Akhiezer function is a solution of 
the formal adjoint equation 

^+(91,. -L,) = (3.13) 

As in the case of differential operators, here and below the right action of a difference 
operator is defined as formal adjoint action, i.e., /^T = r^^/+. 

4 Integrable hierarchies 

In its original form equations (I2.10p . (|2.17p is just an infinite system of partial differ- 
ential equation for an infinite number of coefficients of all the operators, depending 
on infinite number of independent variables called "times" . Of course, restricting to 
a finite number of variables one gets an equation or a finite number of equations for a 
finite number of variables. Some of them are fundamental equations of mathematical 
physics, and as such deserve special interest. That is true for all three basic equations 
mentioned above, that is KP, 2D Toda and BDHE. Our next goal is to present the 
hierarchies of these equations in the form of commuting fiows on a certain "phase 
spaces" that are spaces of pseudodifferential or pseudodifference operators. This form 
is due to Sato and his coauthors [llj . 

KP hierarchy 

Let O be a linear space of a formal pseudodifferential operators in the variable 2;, 
i.e., formal series 

00 

^^(^)^-' (4.1) 

s=-N 

By definition the coefficient vi at 9^7^ in (|4.ip is called the residue of T> 

vi-.^vesgV. (4.2) 
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The commutator relations dx ■ v{x) = Vx{x) + v{x)dx and • v{x) = v{x)dx^ — 
Vx{x)dx^ + Vxx{x)dx^ define on O a structure of associative ring. For any pseudo- 
differential operator T> its differential part is defined as the unique differential oper- 
ator such that 2? — 1)+ = V- = 0{d~^), i.e., for V as in (|4.ip its differential part is 
equal to 



V+= v,{x)dx'^ (4.3) 

s=-N 

The KP hierarchy is defined on the space V of monic pseudodifferential operators of 
order 1, i.e., of the operators of the form 

oo 

C = dx + J2vs{x)d-^ (4.4) 

s=l 



Proposition 4.5 The equations 

d,C=[CX,C] (4.6) 
define commuting flows on the space V . 

Proof. The left-hand side of equation (|4.6p is a pseudodifferential operator diC = 
^g:^i{diVs)d^'' of order at most —1. Therefore, (|4.6p is well-defined if and only if 
the right-hand side is a pseudodifferential operator of order at most —1. To show 
this, notice, that the identity = implies = —[D_,C]. Be definition 

£!_ is an operator of order at most —1. Hence, [ilL,/^] is also of order at most —1. 

For the proof of the second statement of the proposition it is necessary to show 
that equations (|4.6|) imply the equation 

[d, - C\, dj - C?+] = d,C'+ - d,C\ + = (4.7) 

The left-hand side of (|4.7|) is a differential operator. Therefore, in order to show 
that it vanish, it is enough to show that it is a pseudodifferential operator of order 
at most —1. From (j4.6|) it follows that di£^ = [C\,C^] Then using the the identity 
[C\a] = we have 

= [/:V, - d,cl = [c\cl] + o{d-^) = [c^+, c_] + o{d-^) (4.8) 

Similarly, 

[C^,CP^] ^ [C\,C_] [C\^C_]+0{d-^) (4.9) 

Substituting (|4.8p , (|4.9p into (|4.7p completes the proof of the proposition. 

The operator C\ has the form dx — u{x,y), with u — ~1v\ where v\ is the 
coefficient at dx^ of £, i.e., v\ — resg L. Equations (|4.7p with j = 2 have the form 

dt^u = {dy -dl^u, Cl\ ^-{dy-dl^ u, C"J] = 2dxF^ , (4.10) 

where 

F„ :=res9/:™. 
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Important remark At first glance the system (|4.10l) looks like a system of com- 
muting evolution equations, but it is not. The right-hand side of (I4.10p are universal 
differential polynomials in Vi. In general there is no way to reconstruct from one 
function u{x,y) an infinite set of functions Vi{x) of one variable. It can be done 
only under ceratin assumptions. In [41] that was done in the case when u{x, y) is a 
periodic function of the variables x and y. To some extend the main part in the proof 
of the first case of Welter's conjecture can be seen as the proof of the equivalence of 
(|4.6|) and (|4.10l) in the case when u is as in the statement of Theorem 11.61 

For further use let us present some other basic notations and construction. The 
first one is the notion of wave function. 

Lemma 4.11 Let C he a monic pseudodifferential operator of the form ^4-4\ l- Then 
the equation Ctl) = kip has a unique solution of the form 



i^ = e'^(l + f2U^)k-^) (4.12) 

^ s=l 



normalized by the condition S,s{0) = 0. 

The proof is elementary. Substituting (|4.12p into the equation gives a system of 
equations having the form px£,s — Rsivk^^i) with fc, s' < s. Therefore, they uniquely 
define ,if the initial conditions are fixed. 

The wave function is then define the wave operator 

oo 

$ = l + ^^,(x)5r (4.13) 

s=l 

by the equation i/j — $e'^'^. Notice, that the last equation implies 

£ = $ . • (4.14) 
The formal dual wave function is given by the formula 

1/^+ = 6-'=^=^ (l + £ ax)k-") ■■= e-'=^$-i (4.15) 

^ s=l ^ 

is a solution of the formal adjoint equation tp^ C = kip^ 

The defining property of the dual wave function are equations that we proved for 
the dual Baker- Akhiezer function in the previous section. Namely, 

Lemma 4.16 Let tp he a wave function and its dual. Then the equations 

resk{i^+{d^i^))dk = 0, n = 0,l,... (4.17) 

hold. 

The proof is a direct corollary of the identity 

resfc (e-^-"2?i) (Pae'^") dk = resg {V2V1) , (4.18) 

which holds for any pair of pseudodifferential operators (for details see pTl [15] V 

In the same way one can show that the product of the wave function and its dual 
is a generating series for the right-hand sides of the hierarchy (|4.10p . 
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Lemma 4.19 The coefficients of the expansion 

oo 

iP+tP^l + J2Jsk~" (4.20) 

s=2 

are given by Jn+i ~ Fn = resg C" . 

Proof. From the definition of C it follows that 

resfc dk = resfc dk = J„+i. (4.21) 

On the other hand, using the identity (|4.18l) we get 

reSfc(V'+i;"V') dk = reSfc (e^'^^^-i) (i:"$e'=^) dk = lesg = F^. (4.22) 
The lemma is proved. 

2D Toda hierarchy 

In the two-point case there are two sets of continuous variables and one discrete 
variable which we denote by x. It is instructive enough to consider the hierarchy of 
equations corresponding to one set of continuous times associated with one marked 
point. In this subsection we present the definition of the hierarchy of the differential- 
difference equations (|2.17p in the form of the commuting fiows on the space V of the 
pseudodifference operators of the form 

oo 

C = T + J2^s{x)T-', T = e^== (4.23) 

s=0 

In the ring of the pseudodifference operators 

oo 

J2 M^)T~' (4.24) 

s=-N 

the notion of the residue as follows: 

resrV-.^vo (4.25) 

For any pseudodifferential operator V its positive part is defined as the difference 
operator such that V- ■.^V-V+ = 0{T-^, i.e., if V is as in (|424|) . then 

-1 

V+ := Ws{x)T-^ (4.26) 

s=-N 

Proposition 4.27 The equations 

d,C^[C\,C] (4.28) 
define commuting flows on the space V . 



22 



The proof of the first statement goes along the same hnes as in the case of KP 
hierarchy. The proof of the second statement that (I4.28[) implies 



[d,~ C\,d,- C\]^Q (4.29) 

is also identical. The first operator £+ is of the form C+ ~T — u with u = wq. The 
equation (|4.28p for i = 1 gives dtu = —wi, where wi = resT CT. Here and below 
t = ti. For further use, let us present the equation 

dtF^ - (1 - T)F^, (4.30) 

where 

Fra=TesTC"', ^resT-C^T, 

which directly follows from the comparison of residues of two side of the equality 
dtC™ = [C+,Cm]. The commutativity equations (|4.29p imply that the evolution of 
u with respect to all the other times 

dt„u = -iT-l)F^,^ = -dtFm (4.31) 

As in the KP case, in general the last equations can not be regarded as well-defined 
hierarchy on the space of one function u{x,t) because the definition of F^^ involves 
other coefficients of C The main part of the proof of the second case of Welters' 
conjecture can be seen as a reconstruction of C in terms of u under the assumption 
of Theorem [Til 

We conclude this section by providing a necessary definitions and identities, which 
are just discrete analog of that above. Namely, the wave function is a solution of the 
equation Cip — kip of the form 

i;^k^fl + Y,isix)k-') (4.32) 

It defines a unique wave operator by tlie equation 

oo 

tlj = ^k\ $ = l + ^V?,(x)T-^ (4.33) 

s=l 

Then, the dual wave function is defined by the left action of the operator 

i/j^ = Recall that the left action of a pseudodifference operator is the formal 

adjoint action under which the left action of T on a function / is (fT) — T~^f. 

Lemma 4.34 The coefficient of the product 

oo 

^+^=.l + J2Js{Z,t)k-' (4.35) 

s=l 

are equal to Jn = Fn = res^ ■ 

Proof. From the definition of L it follows that 

resfc (V'+(/:"V')) k^^dk = reSfc (V'+fc"^') k'^dk = J„. (4.36) 
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On the other hand, using the identity 

resfe (k-^'Vi) (Pafc^) k'^dk = resy (V2V1) (4.37) 
which is the 2D Toda analogue of (|4.18p . we get 

resk{^p+C"i^)k-^dk = resfe (fc-^^-^) k'^dk = resy >C" = F„. (4.38) 

Therefore, = J„ and the lemma is proved. 

5 Commuting differential and difference operators. 

In the previous section hierarchies of the KP and 2D Toda equations were defined as 
systems of commuting flows on the spaces of pseudodifferential or pseudodifference 
operators, respectively. Consider now the subspace On C of operators whose n-th 
power is a differential (difference) operator L„, i.e., £" = i„ or equivalently — 0. 
The latter directly implies that dt^C ~ 0. In other words the subspace On is the 
subspace of stationary points of the n-th flow of the hierarchy. It has finite functional 
dimension and can be simply identified with the space of all monic differential (dif- 
ference) operators because any such operator L„ uniquely defines the corresponding 

1 /n 

pseudodifferential C = Ln ■ The subspace On is invariant with respect to all the 
other flows. Their restriction on On is a closed system of evolution equations on 
a space of finite-number of unknown functions and can be represented in the form 
diLn = [L^''",L„]. For n = 2 the corresponding reduction of the KP hierarchy is 
equivalent to the hierarchy of the KdV equation 4ut = dtuux + Uxxx- An attempt to 
find explicit periodic solutions of the KdV equation had led Novikov in to the idea to 
consider further reduction to stationary points of one of the "higher" KdV fiows. In 
terms of the original KP hierarchy that is a subspace stationary for two flows of the 
hierarchy (or two linear combinations of basic flows). The corresponding subspace 
is the space of differential order n monic ordinary differential operator L„ such that 
there exists operator Lm commuting with L„ of order m (not multiple of n), i.e., 
the space of solutions of a system (|1.3ip . As it was mentioned in the introduction, 
the problem of classification of commuting ordinary differential operators as pure 
algebraic problem was consider in remarkable works by Burchnall and Chaundy [8] . 

Briefly the key points of their proof of the statement that a pair of such operators 
is always satisfy algebraic relation 

i?(L„,L,„) = 0. (5.1) 

are the following. The commutativity of L„ and L„i implies that the space V^(A) of 
solutions of the ordinary linear equation L„j/(x) — Xy{x) is invariant with respect to 
the operator Lm- The matrix elements of the corresponding finite dimensional 
linear operator Ljn{X) 

Lm\v{X) = im(A) : V{X) ^ V{X) (5.2) 

in the canonical basis Ci{x,X,xo) G C{X), Ci{x, X,xo)\x=xo = ^ij^ are polynomial 
functions in the variable A. They depend on the choice of the normalization point 
X — xq, i.e., Lll^ — L'^-i^lX, xq) . The characteristic polynomial 

R{X,^l)^det{^l-L%{X,xo)) (5.3) 
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is a polynomial in both variables A and and does not depend on xq. 
According to the property of characteristic polynomials we have 

R(Ln,Lm)y{x,\) = 0. 

Notice, that R{Ln,Lm) is an ordinary differential operator. Therefore, if it is not 
equal to zero then its kernel is finite dimensional. Hence, the last equation valid for 
all A implies (|5.ip . and the first statement of [S] is proved. 

The equation i?(A, /i) = defines affine part of an algebraic curve. Let us 
show that it is always compactified by one smooth point Pq. Indeed the equation 
Lni' = k^ij) has always a unique formal wave solution, i.e., a solution of the form 
(|4.12l) normalized by the conditions ^s(0)=0. Moreover, any solution of the latter 
equation of the form e^^ ■ (Laurent series in fc~^) is equal to V'(2;, k)c{k), where c(fc) 
is a constant Laurent series. The operator Lm commutes with L„, therefore Lmtp is 
also a solution to the same equation. Hence, there exists a Laurent series 

oo 

a™(fc) = fc™ + am^sk-' (5.4) 

s= — m+l 

such that Cmip = a{k)ip{x, k), i.e., V' is a formal common eigenfunction of the oper- 
ators Ln, Ljn- That implies the following expansion of the characteristic equation at 
infinity A — > infty: 

ri-l 

R{X,t^) = l[{f^-a{k,)), k^ = X. (5.5) 

4=0 

Now we are ready to explain a role of the condition under which Burchnall and 
Chaundy where able to make the next step. Namely, the condition that orders of 
operators are co-prime. The leading coefficient of a{k) is fc™. Hence, if {n,m) = 1 
then in the neighborhood of the infinite (and, therefore, almost everywhere else) the 
operator Ln{X) has n-distinct eigenvalues, and is diagonalizable, i.e., for each generic 
point P = (A, /i) G F there is a unique eigenfunction ipix, P',xo) of the operators 
Ln,Lm normalized by the condition ^{xq, P;xo) = 1. It can be written as 

n-l 

ip{x,P;xo) = 'Yhi{P,xo)ci{x,X;xo), hQ{P,XQ) = l, (5.6) 

4 = 

where are canonical basis of solution to the equation L„y = Xy defined above 
and hi are coordinates of the eigenvector of the matrix CmiX). They are rational 
expressions in A and /i, and, therefore are meromorphic functions of P G F (if F 
is smooth, otherwise they become meromorphic on an normalization of F). The 
functions Ci, as solutions of the initial value problem, are entire function of the 
variable A. Hence, ip in an affine part of F is a meromorphic function with poles 
that are independent of x (but depend on the normalization point x = xq). If F is 
smooth than their number is equal to the genus g of F. By definition of the canonical 
basis we have that ipx{x, P)i/j~^{x, P)\x=xo h-iiP^ ^o)- The asymptotic of hi can be 
easy found using the formal wave solution. It equals hi ~ k + (0(fc~^)). Therefore 

tp = exp ^ J^^ hi{x, P)dx^ has at Pq exponential singularity and is a Baker- Akhiezer 

function (with the shift of x by xq). 
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Theorem 5.7 |3 [Mil There is a natural correspondence 



(5.8) 



between regular at x = commutative rings A of ordinary linear differential opera- 
tors containing a pair of monic operators of co-prime orders, and sets of algebraic- 
geometrical data {r,Po, where T is an algebraic curve with a fixed first 
jet [k~^]i of a local coordinate in the neighborhood of a smooth point Pq G F and 
J- is a torsion-free rank 1 sheaf on T such that 



The correspondence becomes one-to-one if the rings A are considered modulo conju- 
gation A! — g{x)Ag~'^(x). 

Note that in [29l |30l [8] the main attention was paid to the generic case of the com- 
mutative rings corresponding to smooth algebraic curves. The invariant formulation 
of the correspondence given above is due to Mumford [48] . 

The algebraic curve T is called the spectral curve of A. The ring A is isomorphic 
to the ring A(r,Po) of meromorphic functions on F with the only pole at the point 
Pq. The isomorphism is defined by the equation 



Here "00 is a common eigenfunction of the commuting operators. At x = it is a 
section of the sheaf J- ® 0{—Pq). 

Remark. As we have seen above, the construction of the correspondence (|5.8|) 
depends on a choice of initial point xo = 0. The spectral curve and the sheaf F are 
defined by the evaluations of the coefficients of generators of A and a finite number 
of their derivatives at the initial point. In fact, the spectral curve is independent on 
the choice of xq, but the sheaf does depend on it, i.e., J- — T^^- 

Using the shift of the initial point it is easy to show that the correspondence (|5.8p 
extends to the commutative rings of operators whose coefficients are meromorphic 
functions of x at x = 0. The rings of operators having poles at x = correspond to 
sheaves for which the condition (|5.9|) is violated. 

Remark. In their original paper Burchnall and Chaundy stressed that there is no 
approach to a classification of commutative differential operators whose ordered are 
not co-prime. The classification of commutative rings of ordinary differential opera- 
tors was completed in [32, , where it was shown that a maximal ring A of commuting 
differential operators is uniquely defined by an algebraic curve with marked point, 
the first jet of local coordinate at the marked point, and if the curve is smooth by 
the rank k and degree rg vector bundle. In addition it depends on r — 1 arbitrary 
functions of one variable. Here k is the rank of A defined as the greatest common 
divisor of the orders of commuting operators. 

Commuting difference operators 

A theory of commuting difference operators containing a pair of operators of co-prime 
orders was developed in [48l [31] . It is analogous to the theory of rank 1 commuting 
(Relatively recently this theory was generalized to the case of commuting difference 



ij°(r,j-) = iji(r,J') = 0. 



(5.9) 



LaV'o = aV'o, La^A, ae A(r, Pq). 



(5.10) 
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operators of arbitrary rank in [10] ■) For further use we present here the classification 
of commutative differential operators of the form 

ri-l 

L„ = r" + ^u,(x)r (5.fi) 

s=l 

Theorem 5.12 (\4S[ \31^ ) Let A be a maximum commutative ring of ordinary differ- 
ence operators of the form h5.ll]) containing a pair of operators of co-prime orderes. 
Then there is an irreducible algebraic curve T, such that the ring is isomorphic 
to the ring A(r,P+,P_) of the meromorphic functions on T with the only pole at a 
smooth point P_|_, vanishing at another smooth point P_. The ring is uniquely defined 
by a torsion-free rank 1 sheaves T on T such that 

h°{T,T{nP+ -nP^)) = h^{T,T{nP+ - nP^)) = 0. (5.13) 

The correspondence becomes one-to-one if the rings A are considered modulo conju- 
gation A! — g{x)Ag'''^(x). 

Remark. As in the continuous case the construction of the correspondence depends 
on a choice of initial point xo = Q. The spectral curve and the sheaf T are defined 
by the evaluations of the coefficients of generators of ^ at a finite number of points 
of the form xo + n. In fact, the spectral curve is independent on the choice of xq, 
but the sheaf does depend on it, i.e., T = J-xq- 

Using the shift of the initial point it is easy to show that the correspondence (|5.8p 
extends to the commutative rings of operators whose coefficients are meromorphic 
functions of x. The rings of operators having poles at a; = correspond to sheaves 
for which the condition (|5.13p for n = is violated. 

6 Proof of Welters' conjecture 

As it was mentioned in the introduction the proof of all the particular cases of Welters' 
trisecant conjecture uses different hierarchies: the KP, the 2D Toda, and BDHE. In 
each case there are some specific difficulties but the main ideas and structures of the 
proof are the same. In all the cases the first step is to construct the wave solution. 
It is necessary to emphasize that it is not a wave solution to the ordinary pseudo- 
differential or pseudodifference operators discussed in Section 4. The corresponding 
wave solutions are defined as formal solutions to a partial differential equation. In 
this case there is no way to define such a solution in a unique way without additional 
assumption on a global structure of the coefficients of the equation. As an instructive 
example we present in this section the proof of the first particular case of Welters' 
conjecture, namely, the proof of Theorem ll.61 

First, we prove the implication {A) — (C). Let t{x^ y) be a holomorphic function 
of the variable x in some open domain Z) G C smoothly depending on a parameter 
y. Suppose that for each y the zeros of r are simple, 

r(x,(2/),y) =0, T,(a;,(y), y) ^ 0. (6.1) 

Lemma 6.2 (\^) If equation (|j.7| ) with the potential u — — 29^ lnr(x, y) has a mero- 
morphic in D solution ip(){x,y), then equations hi. 10]) hold. 
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Proof. Consider the Laurent expansions of ■00 and u in the neighborhood of one of 
the zeros Xi of r: 

o 

+ Vi + Wi{x - Xi) + . . . ; 



[x^ x^Y ^g_3^ 
V'o = - — ^ + Pi+ li{x - Xi) + 5i{x - Xi)'^ + — 

(Ah coefficients in these expansions are smooth functions of the variable y). Substi- 
tution of (|6.3p in (|1.7p gives a system of equations. The first three of them are 

aiXi + 2/3i = 0; + aiUi + = 0; + - ^iXi + aiWi = 0. (6.4) 

Taking the y-derivative of the first equation and using two others we get (|1.10l) . 

Let us show that equations (jl.lOp are sufficient for the existence of meromorphic 
wave solutions, i.e., solutions of the form (|1.17p . 

Lemma 6.5 Suppose that equations M.IO]) for the zeros of T{x^y) hold. Then there 
exist meromorphic wave solutions of equation |j. 7| ) that have simple poles at Xi and 
are holomorphic everywhere else. 

Proof. Substitution of ()1.17|) into (|1.7|) gives a recurrent system of equations 

2C+i^dyCs+uis-e: (6.6) 

We are going to prove by induction that this system has meromorphic solutions with 
simple poles at all the zeros Xi of r. 
Let us expand at Xi'. 

= ^ _ ^ + rso + rsi{x - Xi) , (6.7) 

where for brevity we omit the index i in the notations for the coefficients of this ex- 
pansion. Suppose that are defined and equation (j6.6|) has a meromorphic solution. 
Then the right-hand side of (|6.6p has the zero residue a.t x = Xi, i.e., 

lesx, (dy^s + u£,s - = + Virs + 2rsi = (6.8) 

We need to show that the residue of the next equation vanishes also. From (|6.6p it 
follows that the coefficients of the Laurent expansion for Cs+i are equal to 

rs+i = -x^Ts - 2rso, (6.9) 

2rs+i,i = r-so - rsi + Wirg + Vir^o . (6.10) 

These equations imply 

f-s+i + ViTs+i + 2rs+i,i = -rs{xi - 2wi) - Xi{fs - w,rsS + 2rsi) = 0, (6.11) 
and the lemma is proved. 



28 



A-periodic wave solutions 

Our next goal is to fix a translation-invariant normalization of which defines wave 
functions uniquely up to a x-independent factor. It is instructive to consider first the 
case of the periodic potentials u{x + l,y) — u{x,y) (see details in [41]). 
Equations (j6.6|) are solved recursively by the formulae 

e.+i(x,2/) =c,+i(y)+e°+i(x,2/), (6.12) 

y)^l f\Us - e + <s) dx , (6.13) 

where Cs{y) are arbitrary functions of the variable y. Let us show that the period- 
icity condition ^^(x + l,y) = £_s{x,y) defines the functions Cs{y) uniquely up to an 
additive constant. Assume that S,s-i is known and satisfies the condition that the 
corresponding function is periodic. The choice of the function Cs{y) does not affect 
the periodicity property of , but it does affect the periodicity in x of the function 
y)- order to make ^'^^i{x, y) periodic, the function Cs(y) should satisfy the 
linear differential equation 

dyCsiy) + B{y) cs{y) + / y) + u{x, y) y)) dx , (6.14) 



where B{y) — udx. This defines Cs uniquely up to a constant. 

In the general case, when u is quasi-periodic, the normalization of the wave func- 
tions is defined along the same lines. 

Let Yjj = (CU) be the Zariski closure of the group CU = {Ux \ x ^ C} in X. 
Shifting Yu if needed, we may assume, without loss of generality, that Yu is not in 
the singular locus, Yjj (f_ S. Then, for a sufficiently small y, we have Yjj -I- Vy ^ E as 
well. Consider the restriction of the theta-function onto the affine subspace -t- Vy^ 
where :— (the identity component of 7r^^(F[/)), and vr: — > X = / K is the 
universal covering map of X : 

T{z,y)^e{z + Vy), zeC^. (6.15) 

The function u{z, y) — —2di Inr is periodic with respect to the lattice Ajj = A n 
and, for fixed y, has a double pole along the divisor 9 ^(y) — (6 — Vy) fl C''. 

Lemma 6.16 Let equations il.lU\) for zeros of t{Ux-\-z, y) hold and let X be a vector 
of the sublattice Ajy = A n C . Then: 

(i) equation ( [j.7| ) with the potential u{Ux + z, y) has a wave solution of the form 
ip = gkx+k y(j)(jj X + z,y,k) such that the coefficients {z, y) of the formal series 

cb{z, y, k) = e'y (l + f] U^, y) k-'\ (6.17) 



are X-periodic meromorphic functions of the variable z G with a simple pole at 
the divisor Q^{y), 

^,{z + X,y)^U^.,y)^^^; (6.18) 

(a) cf){z, y, k) is unique up to a factor p{z, k) that is du -invariant and holomorphic 
in z, 

(pi{z,y,k) = (l>{z,y,k)p{z,k), dup^O. (6.19) 
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Proof. The functions (,s{z) are defined recursively by the equations 

2dui.+i = dy^s + iu + b)^, - dUs- (6.20) 

A particular solution of the first equation 29c/^i — u + b is given by the formula 

2^? = -2<9(7lnr + (Z,z) 6, (6.21) 

where (Z, z) is a linear form on C'' given by the scalar product of z with a vector I £ 
such that {I, U) = 1. By definition, the vector A is in Yu. Therefore, (Z, A) ^ 0. The 
periodicity condition for Q defines the constant h 

(Z, A)6 = {2du In t{z + \,v) - 2% In r(z, y)) , (6.22) 

which depends only on a choice of the lattice vector A. A change of the potential by 
an additive constant does not affect the results of the previous lemma. Therefore, 
equations (|1.10p are sufficient for the local solvability of (j6.20p in any domain, where 
r(z + Ux, y) has simple zeros, i.e., outside of the set [y) = (Oi - Vy)C^<C'^, where 
01 = n duQ- This set does not contain a 9;7-invariant line because any such line 
is dense in Yu- Therefore, the sheaf Vo of 9[/ -invariant meromorphic functions on 
C^\Q^{y) with poles along the divisor Q^{y) coincides with the sheaf of holomorphic 
^[/-invariant functions. That implies the vanishing of H^iC'^ \ 6;f^(y), Vq) and the 
existence of global meromorphic solutions of (I6.20p which have a simple pole at 
the divisor [y) (see details in [21 [SB]). If ^'"^ fixed, then the general global 
meromorphic solutions are given by the formula = + c^, where the constant of 
integration Cs{z, y) is a holomorphic ^[/-invariant function of the variable z. 

Let us assume, as in the example above, that a A-periodic solution ^^-i is known 
and that it satisfies the condition that there exists a periodic solution of the next 
equation. Let ^*_,_]^ be a solution of (|6.20p for fixed Then it is easy to see that the 
function 

{z, y) = C+i (^, y) + cs{z, y) e?(z, y) + ^9,c,(z, y), (6.23) 

is a solution of (I6.20p for = + c^. A choice of a A-periodic ^[/-invariant function 
Cs(2, y) does not affect the periodicity property of ^s, but it does affect the periodicity 
of the function S^^^j^i- In order to make ^g+i periodic, the function Cs{z,y) should 
satisfy the linear differential equation 

{l,X)dyCs{z,y) = 2C+i(^ + A,y) - 2C+i{z,y) . (6.24) 

This equation, together with an initial condition Cs{z) = Cs(z,0) uniquely defines 
Cs{x, y). The induction step is then completed. We have shown that the ratio of two 
periodic formal series and (j) is y- independent. Therefore, equation (j6.19p . where 
p{z, k) is defined by the evaluation of the both sides at y = 0, holds. The lemma is 
thus proven. 

Corollary 6.25 Let Ai, . . . , A^ be a set of linear independent vectors of the lattice 
Ajj and let zq be a point ofC^. Then, under the assumptions of the previous lemma, 
there is a unique wave solution of equation |_?. 7] ) such that the corresponding formal 
series (j){z , y , k; zq) is quasi-periodic with respect to Kjj , i.e., for A G Ku 

<t>{z + A, y, k; zq) = (j){z, y, k; zq) fix{k) (6.26) 

and satisfies the normalization conditions 

^^x,{k) = l, (/)(zo,0,fc;zo) = 1. (6.27) 
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The proof is identical to that of the part (b) of the Lemma 12 in [SS]- Let us briefly 
present its main steps. As shown above, there exist wave solutions corresponding to 
(f) which are Ai-periodic. Moreover, from the statement (ii) above it follows that for 
any A' G Ac/ 

(j){z + A, y, k) = (f){z, y, k) px{z, k) , (6.28) 

where the coefficients of p\ are ^[/-invariant holomorphic functions. Then the same 
arguments as in [58] show that there exists a 9(7-invariant series f{z, k) with holo- 
morphic in z coefficients and formal series ^\{k) with constant coefficients such that 
the equation 

f{z + A, k)px{z, k) = f{z, k) px{k) (6.29) 

holds. The ambiguity in the choice of / and p corresponds to the multiplication by 
the exponent of a linear form in z vanishing on [/, i.e., 

/'(z,fc)=/(z,fc)e(''W'^), MA(fc)=MA(fc)e(''W'^\ (6(fc),[/) = 0, (6.30) 

where h{k) ~ b^k^^ is a formal series with vector-coefficients that are orthogonal 
to U. The vector U is in general position with respect to the lattice. Therefore, 
the ambiguity can be uniquely fixed by imposing {d — 1) normalizing conditions 
A*Ai(fe) = 1, i > 1 (recall that p,\^{k) = 1 by construction). 

The formal series /</) is quasi-periodic and its multipliers satisfy (|6.27l) . Then, by 
that properties it is defined uniquely up to a factor which is constant in z and y. 
Therefore, for the unique definition of 4>q it is enough to fix its evaluation at zq and 
y = 0. The corollary is proved. 

The spectral curve 

The next goal is to show that A-periodic wave solutions of equation (jl.7p . with u as 
in (jl.8p . are common eigenfunctions of rings of commuting operators. 

Note that a simple shift z — z -f- Z, where Z ^ E, gives A-periodic wave solutions 
with meromorphic coefficients along the affine subspaces Z -I- C^. Theses A-periodic 
wave solutions are related to each other by ^[/-invariant factor. Therefore choosing, 
in the neighborhood of any Z ^ E, a hyperplane orthogonal to the vector U and 
fixing initial data on this hyperplane at y = 0, we define the corresponding series 
(f){z + Z, y, k) as a local meromorphic function of Z and the global meromorphic 
function of z. 

Lemma 6.31 Let the assumptions of Theorem \1.6\ hold. Then there is a unique 
pseudodifferential operator 

oo 

CiZ, d,) =d, + Y. ^siZ)d-' (6.32) 

such that 

C{Ux + Vy + Z,d^)iJj^k^, (6.33) 

where ip = e^^^^ ^(piUx + Z,y,k) is a \-periodic solution of 7| ). The coefficients 
Ws{Z) of C are meromorphic functions on the abelian variety X with poles along the 
divisor Q. 
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Proof. Let be a A-periodic wave solution. The substitution of (|6.17l) in (|6.33p gives 
a system of equations that recursively define Ws{Z,y) as differential polynomials 
in ^siZ,y). The coefhcients of ip are local meromorphic functions of Z, but the 
coefficients of C are well-defined global meromorphic functions of on \ S, because 
different A-periodic wave solutions are related to each other by Scz-invariant factor, 
which does not affect C. The singular locus is of codimension > 2. Then Hartogs' 
holomorphic extension theorem implies that Ws{Z,y) can be extended to a global 
meromorphic function on C^. 

The translational invariance of u implies the translational invariance of the A- 
periodic wave solutions. Indeed, for any constant s the series (j}{V s + Z,y — s,k) 
and 4>{Z, y, k) correspond to A-periodic solutions of the same equation. Therefore, 
they coincide up to a 9t/-invariant factor. This factor does not affect C Hence, 
Ws{Z,y) = Ws{Vy + Z). 

The A-periodic wave functions corresponding to Z and Z + X' for any A' G A are 
also related to each other by a 9c/-invariant factor: 

du {MZ + \\y,k)<p-\Z,y,k)) =0. (6.34) 

Hence, Ws are periodic with respect to A and therefore are meromorphic functions 
on the abelian variety X. The lemma is proved. 

Consider now the differential parts of the pseudodifferential operators C"^. Let 
£!P be the differential operator such that £!" = - = F^d-^ + 0(d-^). The 
leading coefficient Fm of is the residue of £™: 

F„, = resa (6.35) 

From the construction of C it follows that [dy — d'^ + u, C"] = 0. Hence, 

[dy ~dl+ u, = - [dy -dl + u, ri] = 2d^Fra (6.36) 

(compare with ()4.10p ). The functions F„i are differential polynomials in the coeffi- 
cients Ws of C. Hence, Fm{Z) are meromorphic functions on X. Next statement is 
crucial for the proof of the existence of commuting differential operators associated 
with u. 

Lemma 6.37 The abelian functions Fm have at most the second order pole on the 
divisor Q. 

Proof. We need a few more standard constructions from the KP theory. If ■0 is as in 
Lemma 13. 8[ then there exists a unique pseudodifferential operator $ such that 

oo 

^ = $gfcx+fc^y^ $ = l + ^(^,(t/a;-f Z,y)9-^ (6.38) 

s=l 

The coefficients of $ are universal differential polynomials on Therefore, fs{z + 
Z, y) is a global meromorphic function of z G and a local meromorphic function 
of Z ^ S. Note that C = $((9^) 

Consider the dual wave function defined by the left action of the operator 
ij}^ = (e~'^^~'^ ^)$~^. Recall that the left action of a pseudodifferential operator 
is the formal adjoint action under which the left action of dx on a function / is 
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ifdx) = —dxf- If '0 is a formal wave solution of (|1.7p . then is a solution of the 
adjoint equation 

{-dy -dl+ u)iA = 0. (6.39) 

The same arguments, as before, prove that if equations (|1.10p for poles of u hold 
then 5+ have simple poles at the poles of u. Therefore, if "0 is as in Lemma I6.16[ 
then the dual wave solution is of the form ^+ — q-'^^-^ y<j)'^{Ux + ?/, fc), where 
the coefficients £,t{z + Z, y) of the formal series 

0+(z + Z, y, k) = e-'^y (l + £ + y) fc-^) (6.40) 

^ s=l ' 

are A-periodic meromorphic functions of the variable z G with a simple pole at 
the divisor (y). 

The ambiguity in the definition of ip does not affect the product 

= (^e-''^-'''y<S>-^^ (^$e'=^+'='^) . (6.41) 

Therefore, although each factor is only a local meromorphic function on \ S, the 
coefficients Js of the product 

oo 

V'+V'-0+(Z,y,fc)0(Z,y,fc) = l + ^J,(Z,y)fc-^ (6.42) 

are global meromorphic functions of Z. Moreover, the translational invariance of u 
implies that they have the form Js{Z,y) = Js{Z + Vy). Each of the factors in the 
left-hand side of (|6.42p has a simple pole on — Vy. Hence, Js{Z) is a meromorphic 
function on X with a second order pole at 8. According to Lemma [4.19[ we have 
Fn — Jn+1- That completes the proof of the lemma. 

Let F be a linear space generated by {Em, m = 0, 1, . . .}, where we set Fq — 1. It 
is a subspace of the 2^-dimensional space of the abelian functions that have at most 
second order pole at 0. Therefore, for all but g = dimF positive integers n, there 
exist constants Ci,„ such that 

n-l 

F„{Z) + J2c^.nF^{Z)^0. (6.43) 

i=0 

Let / denote the subset of integers n for which there are no such constants. We call 
this subset the gap sequence. 

Lemma 6.44 Let C be the pseudodifferential operator corresponding to a X-periodic 
wave function tjj constructed above. Then, for the differential operators 

Ln = Cl + Y, c^^rXT' nil, (6.45) 

i=0 

the equations 

oo 

L„0 = a„(fc)V, a„(fc) = fc"+^a,,„fc"-^ (6.46) 

s=l 

where Os^n o-fs constants, hold. 
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Proof. First note that from (I6.36P it follows that 



[dy -dl + u, L„] = 0. (6.47) 

Hence, if -0 is a A-periodic wave solution of (|1.7p corresponding to Z ^ Y,, then 
Lnip is also a formal solution of the same equation. That implies the equation 
L^tp = o,n{Z, k)ip, where a is 9^/- invariant. The ambiguity in the definition of il) does 
not affect a„. Therefore, the coefhcients of a„ are well-defined global meromorphic 
functions on \ S. The djj- invariance of a„ implies that a„, as a function of 
is holomorphic outside of the locus. Hence it has an extension to a holomorphic 
function on C^. Equations (|6.34l) imply that a„ is periodic with respect to the lattice 
A. Hence an is Z-independent. Note that Os^n — Cs,n, s < n. The lemma is proved. 

The operator Lm can be regarded as a Z ^ E-parametric family of ordinary 
differential operators whose coefficients have the form 

m 

Li^^^ + Y,^^,^iUx + Z)^:^-\ mil. (6.48) 

i=l 



Corollary 6.49 The operators commute with each other, 

[L,f,i,fj-0, Z^E. (6.50) 

From (|6.46l) it follows that , L^]ijj ~ 0. The commutator is an ordinary differential 
operator. Hence, the last equation implies (|6.50p . 

Lemma 6.51 Let , Z ^ Y,, be a commutative ring of ordinary differential op- 
erators spanned by the operators if. Then there is an irreducible algebraic curve 
F of arithmetic genus g — dimF such that is isomorphic to the ring A(F, Pq) 
of the meromorphic functions on F with the only pole at a .smooth point Pq. The 
correspondence Z — )■ A^ defines a holomorphic imbedding of X into the space of 
torsion-free rank 1 sheaves J- onT 

j:X\T.< — >Plc(F). (6.52) 

Proof. In order to get the statement of the theorem as a direct corollary of Theorem 
5.1, it remains only to show that the ring A^ is maximal. Recall, that a commutative 
ring A of linear ordinary differential operators is called maximal if it is not contained 
in any bigger commutative ring. Let us show that for a generic Z the ring A^ is 
maximal. Suppose that it is not. Then there exits a G I, where / is the gap sequence 
defined above, such that for each Z ^ Y, there exists an operator if of order a 
which commutes with Lf^,n ^ /. Therefore, it commutes with C. A differential 
operator commuting with C up to the order 0(1) can be represented in the form 
La — J2m<a'^i,<^(^)^+^ whcrc Ci^a{Z) are 9i-invariant functions of Z. It commutes 
with C if and only if 

n-l 

Fa{Z) + J2c^AZ)F^{Z)^0, ^uc^,a = 0. (6.53) 

1=0 

Note the difference between (j6.43p and (j6.53|) . In the first equation the coefficients 
Ci^n are constants. The A-periodic wave solution of equation (|1.7|) is a common 
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eigenfunction of all commuting operators, i.e., L^'ip = aa{Z,k)ip, where Ua = k" + 
'^s,a(-^)fc"~* is 9i-invariant. The same arguments as those used in the proof 
of equation (j6.46p show that the eigenvalue Ua is independent. We have Us^a = 
Cs,Q, s < a. Therefore, the coefficients in (|6.53p are Z- independent. That contradicts 
the assumption that a ^ I. The lemma is proved. 

Our next goal is to prove finally the global existence of the wave function. 

Lemma 6.54 Let the assumptions of the Theorem \1.19\ hold. Then there exists a 
common eigenfunction of the corresponding commuting operators Lf^ of the form 
ijj = e^^4){Ux + Z, k) such that the coefficients of the formal series 

oc 

^iZ,k)^l + J2UZ)k-' (6.55) 

s=l 

are global meromorphic functions with a simple pole at Q. 

Proof. It is instructive to consider first the case when the spectral curve T of the rings 
is smooth. Then, as shown in ( '29' ) , the corresponding common eigenfunction 
of the commuting differential operators (the Baker- Akhiezer function), normalized by 
the condition V'oU=o = 1, is of the form (^ |291 150] ) 

^^JiMP) + U^ + mz)^n,n, (6.56) 
6{Ux + Z)e{A{P) + Z) 

(compare with (|2.24p . Here 6{Z) is the Riemann theta-function constructed with 
the help of the matrix of 6-pcriods of normalized holomorphic differentials on F; 
A : r — > J(r) is the Abel- Jacobi map; Q, is the abelian integral corresponding to the 
second kind meromorphic differential dfi with the only pole of the form dk at the 
marked point Pq and 2'Kitj is the vector of its 5-periods. 

Remark. Let us emphasize, that the formula (j6.56p is not the result of solution 
of some differential equations. It is a direct corollary of analytic properties of the 
Baker- Akhiezer function tJjqIx^P) on the spectral curve. 

The last factors in the numerator and the denominator of ()6.56|) are x- independent. 
Therefore, the function 

J 9{A{P) + UX + Z) ,f,(p) 

ipBA = ^ e ^ (6.57) 

e{Ux + Z) 

is also a common eigenfunction of the commuting operators. 
In the neighborhood of Pq the function ipBA has the form 

i'BA^e''^(l + f:^-l^l±^k-\ k^n, (6.58) 

where Ts{Z) are global holomorphic functions. 

According to Lemma 16.511 we have a holomorphic imbedding Z — j{Z) of A" \ S 
into J(r). Consider the formal series ip = j*4'BA- It is globally well-defined out of 
S. If Z ^ 8, then i{Z) ^ G) (which is the divisor on which the condition ()5.9p is 
violated). Hence, the coefficients of V' are regular out of Q. The singular locus is at 
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least of codimension 2. Hence, using once again Hartogs' arguments we can extend 
ip on X. 

If the spectral curve is singular, we can proceed along the same lines using the 
generalization of (|6.57|) given by the theory of Sato r- function ([52])- Namely, a set 
of algebraic-geometrical data (|5.8|) defines the point of the Sato Grassmannian, and 
therefore, the corresponding r-function; T{t;J-). It is a holomorphic function of the 
variables t — (^1,^2, ■ • ^'^d is a section of a holomorphic line bundle on Pic(r). 

The variable x is identified with the first time of the KP-hierarchy, x — ti. 
Therefore, the formula for the Baker-Akhiezer function corresponding to a point of 
the Grassmannian {[52\) implies that the function ipBA given by the formula 

t[x — k, —yk^, — ifc'^, . . . ; J-) , 
= r(x,0,0,.. V) ' ^'-''^ 

is a common eigenfunction of the commuting operators defined by J-. The rest of 
the arguments proving the lemma are the same, as in the smooth case. 

Lemma 6.60 The linear space F generated by the abelian functions {Fq — l,i^m — 
resale™}, is a subspace of the space H generated by Fq and by the abelian functions 
H,=dud,, \n0{Z). 

Proof. Recall that the functions F„ are abelian functions with at most second order 
pole on 8. Hence, a priori g = dimF < 2^. In order to prove the statement of the 
lemma it is enough to show that Fn — djjQn, where Qn is a meromorphic function 
with a pole along 0. Indeed, if Qn exists, then, for any vector A in the period lattice, 
we have Qn{Z + A) = Qn(,Z) + Cn,\- There is no abelian function with a simple 
pole on Q. Hence, there exists a constant g„ and two gi-dimensional vectors In, I'm 
such that Qn = qn + {In, Z) + (Z^, h{Z)), where h{Z) is a vector with the coordinates 

h, = In 61. Therefore, Fn = (L, U) + {I'n, H{Z)). 

Let il}{x^ Z, k) be the formal Baker-Akhiezer function defined in the previous 
lemma. Then the coefficients fs{Z) of the corresponding wave operator $ (16.38^ 
are global meromorphic functions with poles on Q. 

The left and right action of pseudodifFerential operators are formally adjoint, 

i. e., for any two operators the equality (e~'^'^X'i) (2?2e'^^) — e~^^ (T>iV2e'^^) + 

(e-'^^ {Vse'''')) holds. Here V3 is a pseudodifFerential operator whose coefficients 
are differential polynomials in the coefficients of Vi and • Therefore, from (|6.4ip 
it follows that 

00 / QG 

\ ^ 7-1 7 — -S 1,0/ \ ^ . 



= 1 + F,_ik-' ^l + dJ}^ Qsk~' 1 . (6.61) 

The coefficients of the series Q are differential polynomials in the coefficients of 
the wave operator. Therefore, they are global meromorphic functions of Z with poles 
on 0. Lemma is proved. 

The construction of multivariable Baker-Akhiezer functions presented in Section 
2 for smooth curves is a manifestation of general statement valid for singular spectral 
curves: fiows of the KP hierarchy define deformations of the commutative rings A 
of ordinary linear differential operators. The spectral curve is invariant under these 
flows. For a given spectral curve F the orbits of the KP hierarchy are isomorphic 
to the generalized Jacobian J(F) — Pic*'(r), which is the equivalence classes of zero 
degree divisors on the spectral curve (see details in [58 } [29 l ISO j 152]). 
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As shown in Section 4, the evolution of the potential u is described by equation 
(|4.6[) The first two times of the hierarchy are identified with the variables ti = x,t2 = 
y. Equations (|4.6p identify the space Fi generated by the functions duFn with the 
tangent space of the KP orbit at . Then, from Lemma 6.9 it follows that this 
tangent space is a subspace of the tangent space of the abelian variety X. Hence, 
for any Z ^ S, the orbit of the KP flows of the ring is in X, i.e., it defines an 
holomorphic imbedding: 

iz:JiT)' — >X. (6.62) 

From (16.621) it follows that J{T) is compact. 

The generalized Jacobian of an algebraic curve is compact if and only if the curve 
is smooth ( [14J ) . On a smooth algebraic curve a torsion- free rank 1 sheaf is a line 
bundle, i.e., Pic(r) = J(r)- Then (I6.52p implies that iz is an isomorphism. Note that 
for the Jacobians of smooth algebraic curves the bad locus E is empty ([35]), i.e., the 
imbedding j in (|6.52p is defined everywhere on X and is inverse to iz- Theorem 11.61 
is proved. 

7 Characterization of the Prym varieties 

To begin with let us recall the definition of Prym varieties. An involution a : F — > 
F of a smooth algebraic curve F induces an involution a* : J(F) — > J(F) of the 
Jacobian. The kernel of the map 1 + cr* on J(F) is the sum of a lower-dimensional 
abelian variety, called the Prym variety (the connected component of zero in the 
kernel), and a finite group. The Prym variety naturally has a polarization induced 
by the principal polarization on J(F). However, this polarization is not principal, 
and the Prym variety admits a natural principal polarization if and only if a has at 
most two fixed points on F — this is the case we will concentrate on. 

From the point of view of integrable systems, attempts to prove the analog of 
Novikov's conjecture for the case of Prym varieties of algebraic curves with two 
smooth fixed points of involution were made in [611 [59l [7] . In [61] it was shown that 
Novikov-Veselov (NV) equation provides solution of the characterization problem up 
to possible existence of additional irreducible components. In [591 H] the characteri- 
zations of the Prym varieties in terms of BKP and NV equations were proved only 
under certain additional assumptions. Moreover, in [7] an example of a ppav that is 
not a Prym but for which the theta function gives a solution to the BKP equation 
was constructed. Thus for more than 15 years it was widely accepted that Prym 
varieties can not be characterized with the help if integrable systems. 

In [38] the first author proved that Prym varieties of algebraic curves with two 
smooth fixed points of involution are characterized among all ppavs by the property 
of their theta functions providing explicit formulas for solutions of the integrable 2D 
Schrddinger equation, which is one of the auxiliary linear problems for the Novikov- 
Veselov equation. 

Prym varieties possess generalizations of some properties of Jacobians. In [7] 
Beauville and Debarre, and in ^22, Fay showed that the Kummer images of Prym 
varieties admit a 4-dimensional family of quadrisecant planes (as opposed to a 4- 
dimensional family of trisecant lines for Jacobians). Similarly to the case of Jacobians, 
it was then shown by Debarre in jl2j that the existence of a one-dimensional family 
of quadrisecants characterizes Prym varieties among all ppavs. However, Beauville 
and Debarre in [7] constructed a ppav that is not a Prym but such that its Kummer 
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image has a quadrisecant plane. Thus no analog of the trisecant conjecture for Prym 
varieties was conjectured, and the question of characterizing Prym varieties by a 
finite amount of geometric data (i.e., by polynomial equations for theta functions at 
a finite number of points) remained completely open. 

In |24j S. Grushevsky and the first author proved that Prym varieties of un- 
ramified covers are characterized among all ppavs by the property of their Kummer 
images admitting a symmetric pair of quadrisecant 2-planes. That there exists such 
a symmetric pair of quadrisecant planes for the Kummer image of a Prym variety 
can be deduced from the description of the 4-dimensional family of quadrisecants, 
using the natural involution on the Abel-Prym curve. However, the statement that 
a symmetric pair of quadrisecants in fact characterizes Pryms seems completely un- 
expected. 

The geometric characterization of Prym varieties follows from a characterization 
of Prym varieties among all ppavs by some theta-functional equations, which by using 
Riemann's bilinear addition theorem can be shown to be equivalent to the existence of 
a symmetric pair of quadrisecant planes. In order to obtain such a characterization of 
Prym varieties in P^j a new hierarchy of difference equations, starting from a discrete 
version of the Schrodinger equation was introduced, developed, and studied . The 
hierarchy constructed can be thought of as a discrete analog of the Novikov-Veselov 
hierarchy. 

Theorem 7.1 (Main theorem) An indecomposable principally polarized ahelian 
variety (X, 9) £ Ag lies in the closure of the locus Vg of Prym varieties of unramified 
double covers if and only if there exist vectors A, U,V,W E representing distinct 
points in X, none of them points of order two, and constants ci,C2,c^,wi,W2,w^ G C 
such that one of the following equivalent conditions holds: 
{A) The difference 2D Schrodinger equation 



with 



and 



i^n+l^rn+l ^ Un,m{'>Pn+l,m — V'n.m+l) " 1pn,m — 0, (7. 2) 



0{{n + 1)U + mV + :^W + Z) 9{nU + (m + 1)V + i^W + Z) 
''""eiin+lp +{m+l)V + VW + Z)e{nU + mV + VW + Z) ^ ' ' 



0inU + mV + V„„.W + Z) ^1^2 ^3 (c, c^) , (7.4) 

is satisfied for all Z £ X , where 

1 4_('_lVl+™+l 1 n 

V := Vnm 7, : V := I - V, Cnm ■= C3 [C2 ^ j . (7.5) 



(B) The following identity holds: 



.±i~,'A + U + VtW\ fA+U -V±W 



WiW2iciC2r K y ^ j - ^103(^301)=^ K ^ 2 

+ W2C:i{w3.C2r^K ]-K[ — = , 
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where K: 9 z n. (8[e,0](z)) G C^" is a lifting of the Kummer map M.S\) to the 
universal covering of X . 

(C) The two equations ( one for the top choice of signs everywhere, and one for 
the bottom) 

cf^cl e{z + u - V) e{z -u± w) 9{z + v± w) 

+ c'^^cl e{z -u + V) e{z + u± w) e{z ~v± w) 
= cf 2cT2 0[z -u - V) 9{Z + u± w) e{z + v± w) 

+ 9{Z + U + V)e{Z -U ±W)9{Z -V ±W) (7.6) 

are valid on the theta divisor {Z G X : 6{Z) = 0}. 

A purely geometric restatement of part [B) of this result is as follows. 

Corollary 7.7 (Geometric characterization of Pryms) A ppav {X,6) G Ag lies 
in the closure of the locus of Prym varieties of unramified ( etale ) double covers if and 
only there exist four distinct points pi,p2,P3,P4 G X, none of them points of order 
two, such that the following two quadruples of points on the Kummer variety of X : 

{K{pi + £2^2 + £3^3 + £4^4) I Si e {±1}, £2£3£4 +1} 

and 

{K{pi + £2^2 + £3P3 + £4^4) I £i e {±1}, £2£3£4 = -1} 

are linearly dependent. 

Equivalently, this can be stated as saying that [X, 9) lies in the closure of the Prym 
if and only if there exists a pair of symmetric (under the z M> 2pi — z involution) 
quadrisecants of K{X). 

At first glance the structure of the proof is the same as above. It begins with 
a construction of a wave solution of the discrete analog of 2D Schrodinger equation 
(|7.2[) . But in fact, the hierarchy considered involves essentially a pair of functions and 
is thus essentially a matrix hierarchy, unlike the scalar hierarchy arising for the trise- 
cant case. The argument is very delicate, and involves using the pair of quadrisecant 
conditions to recursively construct a pair of auxiliary solutions (essentially corre- 
sponding to the two components of the kernel, only one of which is the Prym). We 
refer the reader to [24]) for details. 

Our goal for this section is to elaborate on the "only if" part of the statement of 
the theorem, because as a byproduct it gives new identities for theta-function which 
are poorly understood an seems require additional attention. 

Four point Baker- Akhiezer function 

Four-point Baker-Akhiezer function depends on three discrete parameters and, as 
was mentioned in Section 2 gives solution to the BDHE equation. For various choice 
of two linear combination of these variables one obtain various linear equation. In 
[55] (see details in [3T) it was shown that the following choice of the "discrete times" 
gives a a construction of algebraic-geometric 2D difference Schrodinger operators. 

Let r be a smooth algebraic curve of genus g. Fix four points ^'i*,^2^ ^ 
and let D = 71 + • • ■ + 7g be a generic effective divisor on F of degree g. By the 
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Riemann-Roch theorem one computes h^{D + n{P^ — P^) + m{P^ — P^)) = 1, 
for any n,m G Z, and for D generic. We denote by 4'n,7n{P), P G T the unique 
section of this bundle. This means that ipn.m is the unique up to a constant factor 
meromorphic function such that (away from the marked points P^) it has poles only 
at 7s, of multiplicity not greater than the multiplicity of 7^ in _D, while at the points 
P^,P^ (resp. P]~,P2~) the function 'ipn,m has poles (resp. zeros) of orders n and m. 

If we fix local coordinates k^^ in the neighborhoods of marked points (it is cus- 
tomary in the subject to think of marked points as punctures, and thus it is common 
to use coordinates such that k at the marked point is infinite rather than zero) , then 
the Laurent series for ipn,m{P), for P eT near a marked point, has the form 



V^n,™ = fc^"(£c,^(n,m)fc-^Y k = k{P), P^P,^, 



(7.8) 



V'«,™ = fc^™^X^(",m)fc-M, fc = fc(P), P-^P2±. (7.9) 
As it was shown in Section 2 the function ipn.m can be expressed as follows; 



2 (p^^ g(l(P) +nU + mV + Z) n,(p)+mn.^p^ 
^„,™(P)-r_ g^X{P) + Z) 



(7.10) 



where for j = 1,2 the differential (& e H°{Kt + P^ + P^) is of the third kind, 
normalized to have residues ^1 at P^^ and with zero integrals over all the a-cycles, 
and ri' is the corresponding abelian integral; we have the following expression r„m 
is some constant, U = A(Pf ) - A{P:^), V = A(P2~) - A{P^), and 

Z = -J2M1s)+^, (7.11) 

s 

where k is the vector of Riemann constants. 

Change of notation Wc use here notation 9 for the Riemann theta-function of F, 
for later use of 6 for the Prym theta function. 

Theorem 7.12 ( j35j ) The Baker- Akhiezer Junction ipn,m given by formula ([7.i0| ) 
satisfies the following difference equation 



Setup for the Prym construction 

We now assume that the curve F is an algebraic curve endowed with an involution a 
without fixed points; then F is a unramified double cover F — > Fq, where Tq — T /a. 
If F is of genus g = 2(7 + 1, then by Riemann-Hurwitz the genus of Fq is g + 1. From 
now on we assume that 5 > and thus 5 > 1. On F one can choose a basis of cycles 
a^, bi with the canonical matrix of intersections • aj — bi ■ bj — 0, Ui ■ bj = 5ij, < 
hi 2g, such that under the involution a we have (T(ao) = oq, 17(60) — ba, o'{aj) — 
CLg+j, cy{bj) = bg^j, I < j < g. If diOi are normalized holomorphic differentials on F 
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dual to this choice of a-cycles, then the differentials duj = diOj — diOg+j , for j = 1 . . . 5 
are odd, i.e., satisfy a*{duk) = —du^, and we call them the normalized holomorphic 
Prym differentials. The matrix of their 6-periods 

likj = <{ du,, 1 < fc,j <5, (7.14) 

is symmetric, has positive definite imaginary part, and defines the Prym variety 

and the corresponding Prym theta function 

0{z) :^0{z,ll), 

for 2: G C^. We assume that the marked points P^, on F are permuted by the 
involution, i.e., P^ = a{P~). For further use let us fix in addition a third pair of 
points P^ , such that also P^ = a{P^). 

The Abel-Jacobi map F ^ J(F) induces the Abel-Prym map A: F — > P{T) 
(this is the composition of the Abel-Jacobi map A: 7 ^ <^(r) with the projection 
J(F) — 7^(F)). There is a choice of the base point involved in defining the Abel- 
Jacobi map, and thus in the Abel-Prym map; let us choose this base point (such a 
choice is unique up to a point of order two in 7^(F)) in such a way that 

A{P) = -A{a{P)). (7.15) 

Admissible divisors 

An effective divisor on F of degree g — \ — 2g^D — + j2g, is called admissible 
if it satisfies 

[D] + [a{D)] ^Kre J(F) (7.16) 

(where Kr is the canonical class of F), and if moreover H^{D + o'(D)) is generated 
by an even holomorphic differential dll, i.e., that 

dr2(7,) = dl7(cr(7,)) = 0, dn^a*{dn). (7.17) 

Algebraically, what we are saying is the following. The divisors D satisfying (|7.16p 
are the preimage of the point Ky under the map 1 -I- cr, and thus are a translate of 
the subgroup Ker{l -I- cr) C J{T) by some vector. As shown by Mumford [4^, this 
kernel has two components — one of them being the Prym, and the other being the 
translate of the Prym variety by the point of order two corresponding to the cover 
F — > Fq as an element in 7ri(Fo). The existence of an even differential as above picks 
out one of the two components, and the other one is obtained by adding A — (t{A) 
to the divisor of such a differential, for some A. statement. 

Proposition 7.18 For a generic vector Z the zero-divisor D of the function 9 {A{P)+ 
Z) on F is of degree 2g and satisfies the constraints |7.i6p and {I.IT^ , i.e., is admis- 
sible. 
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Remark. S. Grushevsky and the first author had been unable to find a complete 
proof of precisely this statement in the literature. However, both Elham Izadi and 
Roy Smith have independently supplied them with simple proofs of this result, based 
on Mumford's description and results on Prym varieties. As pointed out by a referee, 
this result can also be easily obtained by applying Fay's proposition 4.1 in |21) . In 
[23] independent analytic proof was proposed which also can be seen analytic proof 
of some of Mumford's results. 

Note that the function 9{A{P) + Z) is multi-valued on F, but its zero-divisor is 
well-defined. The arguments identical to that in the standard proof of the inversion 
formula (|7.11l) show that the zero divisor D{Z) := 9{A{P) + Z) is of degree g — 1 = 2g. 

Lemma 7.19 For any pair of points conjugate under the involution a there exists 
a unique differential dflj of the third kind (i.e., a dipole differential with simple poles 
at these points and holomorphic elsewhere), such that it has residues =Fl at these 
points, is odd under a, i.e., satisfies dflj — ^a*{dflj), and such that all of its a- 
periods are integral multiples o/tti, i.e., such a differential dfli exists for a unique 
set of numbers Iq, . . . ,lg ^ Z satisfying 

dnj=TTilk, k = 0,...,g. (7.20) 

Indeed, by Riemann's bilinear relations there exists a unique differential dil of the 
third kind with residues as required, and satisfying dQ — for all k. Note, 
however, that then f^^a*{dQ) is not necessarily zero, as the image a{ak) of the 
loop ak, while homologous to ag+k on F, is not necessarily homologic to Og+fc 
(resp. to oq for a{ao)) on f \ {P^}- Thus each integral a*{dfl), being equal 
to 27ri times the winding number of a{ak) around P^ minus that around Pj~ , is 
equal to 2'Kilk for some Ik G We now subtract from dfl the linear combination 
wi {lodujQ +J2k=i^>':i'^^k + dujg+k)) of even abelian differentials to get the desired 
dflj. 

Theorem 7.21 \24f For a generic D = D{Z) and for each set of integers {n,m,r) 
such that 

n + m + r = mod 2 (7.22) 

the space 

H"iD + n{P+ - Pf) + m{P+ - P^) + r{P+ - P^)) 
is one-dimensional. A basis element of this space is given by 

.p^ u 0{A{P)+nU + mV + rW + Z) miP)+miMP)+rn,iP) 
- Ai„^™,. 9{A{P)+Z) 

(7.23) 

where Qj is the abelian integral corresponding to the differential dflj defined by 
lemma l?. 19[ and U, V, W are the vectors of b-periods of these differentials, i.e., 

2TTiUk ^ f d^i, 2'KiVk ^ f dn2, 2mWk ^ f dil^. (7.24) 

Jbk Jbk Jbk 

The proof is identically the same as the proof of ()2.19|) . It is easy to check that 
the right-hand side of (|7.23|) is a single valued function on F having all the desired 
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properties, and thus it gives a section of the desired bundle. Note that the constraint 
l|7.22p is required due to (I7.20|) . and the uniqueness of up to a constant factor, i.e., 
the one-dimensionahty of the above, is a direct corollary of the Riemann-Roch 
theorem. 

Note that bilinear Riemann identities imply 
2U^A{p-)-A{P+), 2V ^ A{P-) ~ A{P+), 2W ^ A{p-) - A{P+). (7.25) 

Let us compare the definition of 'ipn,m defined for any curve F, with that of ipn.m,ri 
which is only defined for a curve with an involution satisfying a number of conditions. 
To make such a comparison, consider the divisor D = D + P^ of degree g = 2g A- 1, 
and let V'n,m be the corresponding Baker-Akhiezer function. 

Corollary 7.26 For the Baker-Akhiezer Junction ipnm corresponding to the divisor 
D = D + P^ we have 

(7.27) 

where v = Vnm is defined in \7.5^ , i.e., is or 1 so that n + m + v is even. 
Corollary 7.28 If n + m is even, then by formulae jl.lO^ , \7.23^ 



e{A{P) +nU + mV + Z) e{A{Po) + Z) 
e{A{P) + Z) 0{A{Po) + nU + mV + Z) 

g«ri+mr2^ (7.29) 



e{A{P) +nU + mV + Z) ejAjPo) + Z) „,^+_^ 



9{A{P) + Z) 0{A{Pa) + nU + mV + Z) 
where ri — Jp^(dfli — dQi), and we recall that Z — A{D) + k, and Z is its image. 

Remark. This equality, valid for any pair of points P, Pq is a nontrivial identity 
between theta functions. The first author's attempts to derive it directly from the 
Schottky-Jung relations have failed so far. 

Notation For brevity throughout the rest of the paper we use the notation: ipn.m ■— 



Lemma 7.30 \24^ The Baker-Akhiezer function ipn,m given by 

0{A{P) + Un + Vm + u^mW + Z) enn,{p)+mn,{P)+u,,^n,(P) 



9{Un + Vni + V^mW + Z) 9{A{P) + Z) ^(2u^^-i)(niMP+)+ma.,{p+)) ' 

(7.31) 

where Vnm = 1 — Vnm o,s in |y.5[), satisfies the equation (77^, i.e., 



V'n+l,m+l ^n,m('0n+l,m V^njm+l) '4^n^m Oi 

with Un^m as in |y.g| ), ( [7.5[ ), where 

ci = ef^^(^3^), C2 = e"^(^3+), C3=e"^(^^) (7.32) 
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Note that the first and the last factors in the denominator of (I7.3ip correspond to a 
special choice of the normalization constants hn.m,u in (|7.23l) : 

i^n^iP^) = {e{Z + W))-\ Un,n = 0, 

(7 33") 

This normalization implies that for even n + m the difference {ipn+i,m+i — i^n,m) 
equals zero at . At the same time as a corollary of the normalization we get that 
(V'n+i.m ~ V'n.m+i) has no pole at P^ . Hence, these two differences have the same 
analytic properties on F and thus are proportional to each other (the relevant is 
one-dimensional by R,iemann-R,och) . The coefficient of proportionality Unm can be 
found by comparing the singularities of the two functions at . 

The second factor in the denominator of the formula ()7.31|) does not affect equa- 
tion (|7.2p . Hence, the lemma proves the "only if" part of the statement (A) of the 
main theorem for the case of smooth curves. It remains valid under degenerations to 
singular curves which are smooth outside of fixed points Qk which are simple double 
points, i.e., to the curves of type {T,a,Qk}- 

Remark. Equation (|7.2I) as a special reduction of (j7.13p was introduced in [IB]. It 
was shown that equation (I7.13|) implies a five-term equation 

^n+l,'m+l ^nm^n+l.'m— 1 ^n.mV'n— l,m+l ^nra'4^n — l,m — l — ^n,m'^n,m (7.34) 

if and only if it is of the form (|7.2p . A reduction of the algebro-geometric construction 
proposed in |35j in the case of algebraic curves with involution having two fixed points 
was found. It was shown that the corresponding Baker- Akhiezer functions do satisfy 
an equation of the form (|7.2|) . Explicit formulae for the coefficients of the equations in 
terms of Riemann theta-functions were obtained. The fact that the Baker- Akhiezer 
functions and the coefficients of the equations can be expressed in terms of Prym 
theta-functions was first obtained in [23] . 

The statement that 'iljn,m satisfy (|7.34p can be proved directly. Indeed all the 
functions involved in the equation are in 

H"{D + {n + l)P+ -{n- l)Pf + (m + 1)^2+ - {m - 1)P^ + v{P+ - P^)) 

By the Riemann-Roch theorem the dimension of the latter space is 4. Hence, any five 
elements of this space are linearly dependent, and it remains to find the coefficients 
of (I7.34[) by a comparison of singular terms at the points Pf*^ , P2 ■ 

Theorem 7.35 f24^ For any four points A, U, V, W on the image T ^ P{r), and 
any Z G 'P(r) the following equation holds: 

e{z + w)x [e{A + u + v + z)9{z - u) e{z - v) 

- 44 e{A + u -V + z) 9{z - u) e{z + v) 

- 44 e{A-u + v + z) 9{z + u) e{z - v) 

+ 44 o{A-u -V + z) 9{z + u) e{z + v)] = 
= e{A + z)x [e{w + u + v + z)e{z - u) 6{z - v) 

- 44 ^{w + u -V + z)9{z -u) e{z + v) 

- 44 ~u + v + z)9{z + u) e{z - v) 
+ 44 ^{w -u -V + z)9{z + u) e{z + v)]. 



(7.36) 
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To the best of the authors' knowledge equation (|7.36p is a new identity for Prym 
theta-functions. For Z such that 9{W + Z) = it is equivalent to equation (|7.6p with 
the minus sign chosen. The second equation of the pair (|7.6p can be obtained from 
(|7.34|) considered for the odd case, i.e., for n + m = 1 mod 2. Using theta functional 
formulas, it can be shown using ()7.34|) that equation (j7.36|) is equivalent to (|7.2|) . 

8 Abelian solutions of the soliton equations 

In [421 143] the authors introduced a notion of abelian solutions of soliton equations 
which provides a unifying framework the elliptic solutions of these equations and and 
algebraic-geometrical solutions of rank 1 expressible in terms of Riemann (or Prym) 
theta-function. A solution u{x, y, t) of the KP equation is called abelian if it is of the 
form 

-2dl\YiT{Ux + z,y,t) , (8.1) 

where x, y, t G C and z G C" are independent variables, U G C, and for all y, 
t the function T{-,y,t) is a holomorphic section of a line bundle C = C{y,t) on an 
abelian variety X = C"/A, i.e., for all A e A it satisfies the monodromy relations 

r(z + A,y,i) =e°^-^+''^T(z,y,t), for some ax e C , bx ^ bx{y,t) e C . (8.2) 

There are two particular cases in which a complete characterization of the abelian 
solutions has been known for years. The first one is the case n = 1 of elliptic solutions 
of the KP equations. The second case in which a complete characterization of abelian 
solutions is known is the case of indecomposable principally polarized abelian variety 
(ppav). The corresponding 0- function is unique up to normalization, so that Ansatz 
(18. ip takes the form u = —2d^\n6{Ux + Z{y,t) + z). Since the flows commute, 
Z{y,t) must be linear in y and t: u = —2d^ \n9{Ux + Vy + Wt + z) . Besides these 
two cases of abelian solutions with known characterization, another may be worth 
mentioning. Let P be a curve, P € F a smooth point, and tt : F Fq a ramified 
covering map such that the curve Fq has arithmetic genus 50 > and P is a branch 
point of the covering. Let J(F) = Pic^iT) be the (generalized) Jacobian of F, let 
Nm: J(F) — > J(Fo) be the reduced norm map as in [SD], and let 

X = ker(Ar77i)° C J(F) 

be the identity component of the kernel of Nm. Suppose X is compact. By assump- 
tion we have 

dim J(r) - dim A dim,/(Fo) = .go > 0, 

so that A is a proper subvariety of J(F), and the polarization on A induced by 
that on J(F) is not principal, and define the KP flows on PicS^^iT) using the data 

(r,P,C)- 

In general, since for any ro G Z>o the space X]r<ro ^cl/9tr is independent of the 
choice of C, for any C G mp \mp and < r < to (so in particular for r = 1), the r-th 
KP orbit of T is contained in F ® X, and so it gives an abelian solution. Let us call 
this the Prym-like case. An important subcase of it is the quasiperiodic solutions of 
Novikov-Veselov (NV) or BKP hierarchies. 

In the Prym-like case, just as in the NV/BKP case we can put singularities to 
F and Fq in such a way that A remains compact, so it is more general than the 
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KP quasiperiodic solutions. Recall that NV or BKP quasiperiodic solutions can be 
obtained from Prym varieties Prym(r, i.) of curves F with involution l having two 
fixed points. The Riemann theta function of J(r) restricted to a suitable translate 
of Pryni(r, l) becomes the square of another holomorphic function, which defines the 
principal polarization on Prym(r, l). The Prym theta function becomes NV or BKP 
tau function, whose square is a special KP tau function with all even times set to 
zero, so any KP time-translate of it 

• gives an abelian solution of the KP hierarchy with n = dim X being one-half 
the genus g(r) of F, and 

• defines twice the principal polarization on X. 

A natural question is whether these conditions characterize the (time-translates of) 
NV or BKP quasiperiodic solutions. 

Hurwitz' formula tells us that in the Prym-like case n — dim(X) > g{r)/2, where 
the equality holds only in the NV/BKP case. At the moment we have no examples 
of abelian solutions with 1 < n < g(r)/2. 

For simplicity we present here a solution to the classification problem of abelian 
solutions of the KP equation obtained in |42) under an additional assumption on the 
density of the orbit CU mod A in X. 

Theorem 8.3 Let u{x,y,t) be an abelian solution of the KP such that the group 
CU mod A is dense in X . Then there exists a unique algebraic curve F with smooth 
marked point P G F, holomorphic imbedding jo ; X — > J{^) ^.i^d, a torsion-free rank 
1 sheaf J- G Pic^^^(F) where g — g{T) is the arithmetic genus ofT, such that setting 
with the notation j{z) — jo{z) (E) J- 

t{Ux + z,y,t) ^ p{z,y,t)?{x,y,t,0, . . . \ T,P,j{z)) (8.4) 

where T(ti, ^2, ^3, • • • | F, P, J^) is the KP t -function corresponding to the data (F, P, J-), 
and p{z,y,t) ^ satisfies the condition djjp = 0. 

Note that if F is smooth then: 

T{xMMr--\T,P,j{z))^e(ux + Y,V,t,+3{z) I i?(r)) e«(^'*-*--) , (8.5) 

where Vi G C", Q is a quadratic form, and B{T) is the period matrix of F. A 
linearization on J(F) of the nonlinear (y, t)-dynamics for T{z,y,t) indicates the pos- 
sibility of the existence of integrable systems on spaces of theta-functions of higher 
level. A CM system is an example of such a system for n = 1. 

Without the density assumption there are examples in which the KP hierarchy 
has basically no control beyond the closure of the orbit, showing the importance of 
the principal polarization in a Novikov-like conjecture in which a minimal number of 
equation is used to study the nature of X. Having this in mind, we may regard prin- 
cipally polarized Prym-Tjurin varieties |28j as a way to study analogues of Novikov's 
conjecture. 
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